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RUNOFF  CURVE  NUMBERS  FOR  RANGELAND  FROM  LANDSAT  DATA 
1 .  Background 

Within  the  framework  of  a  MSc  program  at  the  Department  of  Hydraulics  and 
Catchment  Hydrology  of  the  Agricultural  University  in  Wageningen,  the  author 
conducted  a  study  at  the  Hydrology  Laboratory  in  Beltsville,  Maryland.  The 
Laboratory  is  part  of  the  Agricultural  Research  Service  (ARS)  of  the  United 
States  Department  of  Agriculture.  The  research  program  at  the  Hydrology 
Laboratory  is  directed  at  developing  methodology  for  studying  hydrologic 
problems  that  encompass  large  areas  of  the  country  and  will  have  a  widespread 
and  universal  application.  The  two  approaches  that  are  emphasized  are 
hydrologic  modeling  and  remote  sensing.  Projects  are  carried  out 
cooperatively  with  the  regional  ARS  hydrology  research  centers,  state  and 
other  federal  agencies.  This  work  contributes  to  the  research  on  the 
hydrology  component  of  the  SPUR  model  (Simulation  of  Production  and 
Utilization  of  Rangelands),  a  physically  based,  rangeland  simulation  model 
developed  by  ARS  to  aid  both  resource  managers  and  researchers. 

SPUR  is  composed  of  five  basic  components:  1)  climate;  2)  hydrology;  3) 
plant;  4)  animal  (both  domestic  and  wildlife);  and  5)  economic.  The  model  is 
driven  by  daily  inputs  of  rainfall,  maximum  and  minimum  temperatures,  solar 
radiation,  and  wind  run,  either  obtained  from  actual  weather  records  or 
stochastically  generated.  Although  specifically  intended  to  furnish  the  SPUR 
model  with  an  input  parameter  on  which  the  water  balance,  erosion,  and  water 
quality  aspects  depend,  i.e.,  the  runoff  curve  number,  this  project  also  has 
significant  potential  for  general  use  by  the  Soil  Conservation  Service  (SCS) 


in  the  western  USA. 
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The  term  rangeland  or  range  refers  to  any  kind  of  grazable  land,  from 
desert,  its  most  marginal  form,  to  grazed  forests,  natural  grasslands,  and 
marshes.  In  the  United  States  there  are  331.8  million  hectares  of  range- 
land.  It  is  because  of  this  great  expanse  of  rangeland,  rather  than  its 
productivity  relative  to  forest  and  croplands,  that  rapid,  low-cost  evaluation 
and  management  techniques  are  required. 

The  runoff  curve  number  (CN)  method  employed  as  a  rainfall  runoff 
relationship  in  this  study  is  widely  used  through  the  spectrum  of  hydrologic 
application  by  federal,  state  and  private  hydrologists.  At  some  point  metric 
units  should  be  used,  however,  in  general  the  philosophy  behind  the  expression 
Min  Rome,  do  as  the  Romans  do”  is  practiced,  mainly  with  the  possible 
applications  of  this  research  in  mind. 
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2.  Introduction 

2.1  Outline 

The  main  objective  is  to  develop  a  procedure  for  estimating  runoff  curve 
numbers  (CN)  for  rangeland  from  existing  Landsat  multispectral  scanner  (MSS) 
information.  The  approach  used  in  this  study  is  to  analyse  remote  sensing 
data  for  small  watersheds  and  to  compare  the  results  to  CN  values  determined 
with  conventional  rainfall  (P)  and  runoff  (Q)  records. 

Earlier  research  has  shown  that  CN  can  be  estimated  using  Landsat  imagery 
and  the  SCS  model  (Bondelid,  et  al.,  1980,  Ragan  and  Jackson,  1980;  Slack  and 
Welch,  1979).  The  utilization  of  Landsat  data  in  a  model  is  generally  time 
and  cost  effective  because  the  Landsat  data  are  in  machine  readable  form. 
Therefore,  digital  computers  can  be  used  for  the  required  computations.  In 
contrast  to  prior  research  which  used  remote  sensing  classified  land  cover  and 
conventional  soils  maps  to  determine  CN,  the  intention  of  this  study  is  to 
relate  Landsat  MSS  data  directly  to  the  measured  CN  of  rangeland  watersheds. 
The  reason  for  omitting  soils  data  from  the  final  model  is  that  soils  maps  are 
generally  not  available  for  vast  areas  of  the  world  as  well  as  for  large  parts 
of  the  western  USA.  Figure  1  shows  areas  for  which  SCS  soils  maps  are 
available  as  of  January  1,  1984. 

A  number  of  rangeland  watershed  records  have  been  collected  and  analysed 
to  develop  basic  P-Q  relations  and  to  get  some  insight  into  the  range  of  CN 
values  and  the  hydrological  behavior  of  the  various  rangelands.  Ten  of  these 
watersheds  were  chosen  for  this  study  to  provide  a  sample  of  data  that  was 
representative  of  the  many  cover-climatic  complexes  found  in  rangelands.  For 
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the  selected  watersheds  Landsat  computer-compatible  tapes  were  obtained  and 
the  spectral  response  was  correlated  with  their  CNs. 


< 

Figure  1.  Availability  of  modern  soil  surveys  as  of  January  1,  198*1. 

Figure  2  presents  a  flow  diagram  of  the  procedures  followed  in  this  research. 
In  the  concluding  chapter  an  outline  and  recommendations  are  given  for  the 
additional  research  needed  to  develop  a  complete  methodology 

2.2  Runoff  curve  numbers 


2.2.1  The  SCS  model 
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In  the  early  1950s  the  Soil  Conservation  Service  (SCS)  of  the  USDA 
developed  a  procedure  for  estimating  stormflow  volumes  from  rainfall  events  on 
small  watersheds.  The  SCS  method  was  designed  to  use  total  daily  (24  hour) 
rainfall.  The  relationships  between  storm  rainfall,  watershed  characteristics 
and  stormflow  are  derived  in  the  so  called  NEH-4 ,  i.e.  The  National 
Engineering  Handbook,  Section  4  (USDA -SCS,  1972). 


OUTLINE  OF  THE  RANGELAND  CURVE  NUMBER  STUDY 


Figure  2.  Flow  diagram  of  the  rangeland  curve  number  study 
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Analyses  of  rainfall  and  runoff  records  indicates  that  the  rainfall 
magnitude  must  be  sufficient  to  satisfy  interception,  depression  storage,  and 
the  infiltration  quantity  before  the  start  of  runoff.  This  threshold  is 
called  the  initial  abstraction,  la. 

After  start  of  runoff,  additional  loss  occurs  mainly  in  the  form  of 
infiltration.  The  total  actual  retention  after  runoff  begins  is  given 
the  symbol  F,  which  increases  up  to  some  maximum  retention,  S.  The  ratio  of 
actual  retention  F  to  maximum  retention  S  is  assumed  to  be  equal  to  the  ratio 
of  runoff  to  rainfall  minus  initial  abstraction,  which  can  be  expressed  as 

Q/(P-Ia)  =  F/S  [l] 

For  the  limit,  as  P  -*■  <»,  F  ->•  S  and  the  ratio  F/S  ■+  1.  The  ratio  of 
Q/(P-Ia)  also  approaches  1,  although  it  can  never  actually  reach  1,  but  for 
all  practical  purposes  the  two  ratios  approach  1  as  ».  When  P  =  la,  F  =  0 
and  the  ratio  of  F/S  =  0.  As  P  becomes  greater  than  la  the  ratio  of  F/S  is 
still  near  zero  and  the  Q/(P-Ia)  ratio  is  also  near  zero.  Since  the 
relationship  holds  at  the  two  end  points,  it  is  assumed  to  hold  for  all 
intermediate  points.  After  runoff  begins,  all  rainfall  becomes  runoff  or 
actual  retention.  That  is. 

(P-Ia)  =  F  +  Q  [2] 

Solving  equations  [l]  and  [2]  for  Q  when  P>Ia  yields 

Q  =  (P-Ia)2  /( (P-Ia)  +  S)  [3] 

and  when  P  <  la , 

Q  =  0 

In  order  to  avoid  estimatimg  both  la  and  S,  the  SCS  expressed  la  in  terms 
of  S  by  the  empirical  relationship 

la  =  0.2S  [4] 


thus  simplifying  equation  [3]  to 
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Q  =  (P-0.2S)2  /  (P+0.8S)  ,  P>0.2S  [5: 

Q  =  0  P_<_0.2S 

Since  S  can  theoretically  range  from  zero  to  infinity,  it  was  transformed  to  a 
scale  from  0  to  100.  The  transformation  takes  the  form 

CN  =  1000  /(S+10)  [6] 

where 

CN  =  runoff  curve  number  or  the  hydrologic  soil-cover-complex  number. 

Curve  numbers  are  reflective  of  land  condition,  and  tables  and  charts  of  CN 
defined  on  soil  type,  land  use,  cover,  and  watershed  moisture  status  are  given 
in  numerous  SCS  and  other  agency  documents. 


Figure  3* 


The  curve  number  nomogram. 
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The  watershed  moisture  status  is  one  of  the  most  influential  factors  in 
CN.  NEH-4  gives  in  its  Table  10.1  conversions  through  three  antecedent 
moisture  classes,  AMC  I,  AMC  II  and  AMC  III,  where  condition  II  is  the 
benchmark  situation  upon  which  a  watershed  is  described,  and  conditions  I  and 
III  are  the  "dry"  and  "wet"  situations  respectively.  An  abbreviated  form  of 
these  relationships  is  given  here  as  Table  1 .  NEH-4  also  gives  suggested 
approximate  five-day  antecedent  rainfall  depths  necessary  to  achieve  these 
three  categories. 


Table  1  .  Relationships  between  CN  and  AMC 


I 

Antecedent  Moisture  Class 

II 

III 

100 

100 

100 

78 

90 

96 

63 

80 

91 

51 

70 

85 

40 

60 

78 

31 

50 

70 

22 

40 

60 

15 

30 

50 

Graphically  the  direct  runoff  is  estimated  from  Figure  3  (NEH-4  Figure 
10.1).  Entering  with  the  total  rainfall  and  CN,  the  runoff  amount  is  read. 

For  example,  a  rainfall  total  of  3  inches  and  a  CN  of  80  gives  a  direct  runoff 
(Q)  of  2  inches. 

Table  2  is  an  excerpt  of  Table  9.1  from  the  NEH-4.  Only  the  portion  for 
pasture  or  range,  meadow,  and  woods  are  presented. 


-9- 


Table  2.  Runoff  curve  numbers  for  hydrologic  soil-cover  complexes 


(Antecendent  moisture 

condition  II,  and  la 

=  0.2 

S). 

Land  use 

Treament 

Hydrologic 

Hydrologic 

soil 

group 

or  practice 

condition 

A 

B 

C 

D 

Pasture 

Poor 

68 

79 

86 

89 

or  range 

Fair 

49 

69 

79 

84 

Good 

39 

61 

74 

80 

Contoured 

Poor 

47 

67 

81 

88 

fl 

Fair 

25 

59 

75 

83 

t! 

Good 

6 

35 

70 

79 

Meadow 

Good 

30 

58 

71 

78 

Woods 

Poor 

45 

66 

77 

83 

Fair 

36 

60 

73 

79 

Good 

25 

55 

70 

77 

The  hydrologic  condition  can  be  determined  by  using  vegetative  condition 
and  air-dry  weight  of  the  cover  (NEH-4  Tables  8.1  and  8.2).  The  appropriate 
Hydrologic  Soil  Group  (HSG)  can  be  determined  by  several  procedures.  Table 
7.1  in  NEH-4  presents  a  list  of  over  4,000  soils  giving  the  HSG.  The  SCS 
Soils-5  file,  put  in  an  interactive  soils  information  system  by  the  U.S.  Army 
Corps  of  Engineers  (1983)  also  lists  the  HSG  for  particular  soil.  Brakensiek 
et  al.  (1984)  presented  a  method  to  estimate  the  HSG  from  textural 


composition . 
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2.2.2  The  statistical  approach 

While  values  of  CN  are  usually  estimated  from  tables  based  on  soils  and 
vegetation,  it  is  also  possible  to  determine  an  actual  CN  from  rainfall  and 
runoff  data.  This  is  possible  through  solution  of  equation  5  for  S  to 

S  =  5(P+2Q-(4Q2+5PQ)1/2)  [7] 

thus  defining  an  S  (and  therefore  a  CN  via  equation  6  )  from  any  P-Q  pair. 
Some  division  exists  in  the  literature  on  the  number  and  kind  of  data  required 
to  define  a  consensus  CN  for  a  watershed.  Two  types  of  P-Q  pairs  can  be  used 
for  the  CN  calculation,  the  natural  events  or  the  ordered  data,  where  the  P 
and  Q  sets  both  in  an  ordered  array  are  matched  up,  with  equation  5  merely 
transforming  the  P  vector  to  the  Q  vector.  The  last  method  was  presented  by 
Hjelmfelt  (1980)  and  preserves  the  original  intended  use  of  the  method  to 
predict  a  given  return  period  Q  from  that  same  return  period  P. 

In  this  study  three  methods  of  CN  determination  are  used  in  the  analysis. 

1.  Mean;  curve  numbers  for  all  the  events  are  calculated  and  the  arithmetic 
mean  is  taken. 

2.  Non-linear  regression;  the  optimum  curve,  describing  the  two  parts  of 
equation  [5],  is  fitted  through  the  P-Q  pairs  using  a  non-linear 
regression  algorithm. 

3.  Relative  storm  size;  only  the  large  rainstorms  are  selected  for  this 
analysis,  using  a  selection  criterion  which  is  described  below. 

For  all  three  methods  both  the  natural  and  ordered  data  were  used  as 
input,  thus  producing  a  total  of  six  CNs  per  watershed.  The  computational 
procedures  or  each  method  are  given  in  section  3*1* 

The  third  method  was  developed  because  small  storms  tend  to  bias  the 
relationship  towards  an  overestimated  CN.  Smith  and  Montgomery  (1980)  and 
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Hjelmfelt  (1982)  suggested  that  "small"  events  should  be  eliminated  from  the 
data  set.  However,  on  a  very  porous  forested  watershed  of  CN  =  40,  a  storm  of 
three  inches  would  be  "small"  while  for  a  near  impervious  parking  lot  of  CN  = 
99,  virtually  any  rainstorm  would  be  "large".  Obviously  there  is  a  need  for  a 
relative  storm  size  definition.  In  order  to  obtain  an  unbiased  storm  size, 
equation  [5]  is  easily  put  into  dimensionless  form  by  dividing  by  S,  thus 
yielding 

Q/S  =  (P/S  -  0.2)2/(PS+0.8) ,  (P/S>  0.2)  [8] 

Q/S  =  0  ,  (P/S  <  0.2 

The  S  values  corresponding  to  the  CN  (in  the  three  AMC)  are  found  to  present  a 


close  correlation  (Hawkins,  et  al.,  1984): 

Sz  =  2.281  Sn  19 1 

and 

SIII  =  0,427  SII  Cl0] 

This  leads  us  to  the  following  rainfall  runoff  equations 

AMC  III  :  Q/Sjj  r  (P/SI];  -0.085)2/(P/Si;[+0.342)  ,  (P/SIX  >  0.085)  till 

Q/SIX  =  0  ,  (P/SI3;  <  0.085) 

AMC  I  Q/Sjj  =  P/SI]:  -  0.456)2/(P/SII+1 .84)  ,  (P/SIX  >  0 .456)  Cl2 3 

Q/SIX  =  0  ,  (P/SIX<  0.456) 


The  equations  [l  l]  and  [12]  are  shown  in  Figure  4,  which  illustrates  CNs  in 
dimensionless  form  for  three  AMCs. 

A  broader  interpretation  of  the  AMC  categories  is  as  "error  bands"  or 
envelopes,  indicating  the  experienced  variability  in  rainfall-runoff 
relationships.  Hjelmfelt  et  al.  (1981)  found  that  established  AMC 
relationships  fairly  accurately  described  the  10  percent,  50  percent,  and  90 
percent  cumulative  probability  of  the  runoff  depth  for  a  given  rainfall  via 
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the  CN  method,  corresponding  to  AMC  III,  II,  and  I  respectively.  An 
illustration  is  given  in  Figure  5. 


Figure  4.  Dimensionless  rainfall-runoff  relationships  for  three  antecedent 
moisture  conditions. 


Figure  5.  Probabilities  associated  with  AMC  categories. 
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Thus  equations  [ll]  and  [l2j  have  associated  runoff  probabilities  which  are 
shown  in  Table  3» 

Table  3*  Runoff  probability  for  P/S* 


P/S 

Prob( Q/S  0) 

Comments 

0 

0 

0.085 

0.10 

AMC  III 

0.20 

0.50 

AMC  II 

0.456 

0.90 

AMC  I 

1 

*S  defined  on  AMC  II 

An  extensive  description  of  this  selection  method  is  in  Appendix  1 .  In 
this  paper  the  criterion  P/S  >  0.456,  or  Prob  (Q/S  >  0)  =  0.90  is  selected  to 
present  a  large  event,  where  S  is  defined  on  AMC  II. 

The  elimination  of  small  and  medium  events  (Appendix  1)  from  the  data  set 
is  an  iterative  procedure  in  which  the  average  Sjj  is  calculated,  the  smallest 
event  is  dropped  if  it  does  not  meet  the  selection  criterion,  a  new  Sjj  is 
calculated,  and  so  forth  until  only  large  events  are  left  in  the  set.  From 
the  remaining  data,  the  average  CN  is  calculated. 

The  disadvantage  of  this  P/S  "selection"  method  is  that  it  does  not 
always  produce  a  CN,  particularly  in  arid  regions,  because  all  the  events  are 
classified  as  "small"  and  therefore  left  out.  Since  of  the  two  remaining 
techniques  the  regression  method  is  the  closest  to  the  P/S  "selection"  method, 
especially  when  calculated  from  ordered  P  and  Q  arrays,  this  method  seems  the 
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most  appropriate  for  the  purpose  of  this  study.  Illustration  of  this  decision 
and  examples  of  the  hydrologic  data  manipulation  are  presented  in  section  3-1 • 

2.3  Landsat  images 

The  first  Earth  Resources  Technology  Satellite,  which  is  now  called  Landsat  1 
was  launched  on  July  23,  1972.  Following  Landsat  1  four  more  Landsats  have 
continued  the  program.  The  satellites  circle  the  globe  in  a  circular,  near 
polar  orbit  so  that  the  same  point  on  the  Earth’s  surface  is  viewed  every  18 
days  at  the  same  time  of  day.  Landsats  1-5  all  carried  a  high  resolution 
sensor,  the  multispectral  scanner  (MSS).  The  Landsat  MSS,  operating  in  the 
visible  and  near  infrared  wavelengths,  continually  scans  the  earth  in  a  185-km 
swath  nominally  perpendicular  to  the  satellite  orbital  track.  The  location  of 
visible  and  infrared  radiation  in  the  spectrum  of  electromagnetic  wavelengths 
is  given  in  Figure  6.  Scanning  is  accomplished  in  the  cross-track  direction 
by  an  oscillating  mirror.  The  forward  motion  of  the  satellite  provides  the 
along-track  progression  of  the  scan  lines.  For  each  mirror  sweep,  six 
adjacent  lines  are  scanned  simultaneously  in  each  of  the  four  spectral  bands 
of  the  electromagnetic  spectrum.  These  four  bands  are: 

1.  Band  4,  the  green  band,  which  covers  the  wavelengths  between  500  and  600 
nanometers  (nm  =  10“^m).  This  band  responds  to  reflected  radiation  from 
sediment-laden  and  shallow  waters. 

2.  Band  5,  the  red  band,  between  600  and  700  nm.  This  band  responds  to  man¬ 
made  features  such  as  urban  and  rural  settlements. 

3.  Band  6,  a  near  infrared  (IR)  band,  between  700  and  800  nm.  This  band 
responds  to  vegetation  and  landforms. 
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4.  Band  7,  another  near  IR  band,  between  800  and  1100  nm.  This  band  also 
responds  to  vegetation,  land  forms,  water,  and  provides  the  best  haze 
penetration . 
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Figure  6.  The  electromagnetic  spectrum. 


The  Landsat  data  are  distributed  in  sections  called  "scenes”.  Each  scene 
is  approximately  1 85  km  wide  by  1 85  km  long.  Each  digital  picture  element 
(pixel)  is  57  m  wide  by  79  tn  long.  The  Landsat  computer  compatible  tapes 
(CCT's)  contain  the  measured  radiance  response  values  in  the  four  bands  for 
each  pixel  in  the  scene.  The  components  of  Landsat  scenes  are  illustrated  in 
Figure  7.  The  rows  of  pixels  are  referred  to  as  "scan  lines”  and  the  columns 
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are  referred  to  as  elements.  The  lines  are  numbered  consecutively  from  north 
to  south  and  the  elements  are  numbered  from  west  to  east. 


I*-  ~~i 

METERS 


Figure  7.  Landsat  scene  components. 


Pattern  recognition  to  distinguish  land  cover  with  Landsat  MSS  is  based 
on  the  premise  that  each  land  cover  type  has  a  specific  range  of  spectral 
responses  in  one  or  more  channels.  To  illustrate  the  possibility  of 
discrimination  by  comparison  of  the  responses  in  different  wavelength  bands, 
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the  typical  relative  response  curves  for  vegetation,  soil,  and  water  are  given 
in  Figure  8.  Figure  9  illustrates  the  principle  of  a  MSS  classifier.  The 
categories  from  Figure  8  are  plotted  in  a  measurement  space  whose  axes  are 
their  band  5  and  band  7  responses.  In  reality  there  is  not  such  a  unique 
measurement  pattern  associated  with  each  category.  Rather,  associated  with 
each  category  is  a  probability  distribution,  due  to  natural  random  variations, 
systematic  seasonal  causes,  atmospheric  conditions,  etc.  In  section  3*2  a 
number  of  MSS  band  combinations  are  described  as  well  as  the  cover  materials 
that  they  emphasize. 


Figure  8.  Relative  response  curves 


-18- 


o 

z 

< 

m 

UJ 

GO 

z 

O 

Q. 

CO 

UJ 

cr 


r 

VE 

3ETAT 

• 

ION 

■ 

« 

SOIL 

• 

water 

_ 1 _ 

RESPONSE.  BANO  5 

Figure  9.  The  categories  vegetation,  soil,  and  water,  and  their  responses  in 
bands  5  and  7. 


A  great  deal  of  this  Landsat  analysis  was  performed  using  the  Office  for 
Remote  Sensing  of  Earth  Resources  (ORSER)  of  Pensylvania  State  University 
computer  software  package  described  by  Borden,  et  al.  (1975).  The  ORSER 
software  provides  substantial  capabilities  for  analyzing  and  classifying 
multispectral  scanner  data.  The  ORSER  software  package  is  fully  operational 
on  the  U.S.  Department  of  Agriculture  (USDA)  IBM  370  computer  system  operated 
by  the  Washington  Computer  Center  (WCC) .  From  a  total  of  approximately  40,  the 
ORSER  programs  that  were  used  in  this  analysis  are  summarized  below. 

TPINFO  Outputs  information  about  contents  of  an  original  multispectral  image 
data  tape  or  SUBSET  data  tape.  The  information  of  interest  in  the 
Identification  (ID)  number,  the  table  of  contents,  the  data  of 
exposure,  and  the  sun  elevation  angle. 
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SUBSET  Allow  the  user  to  specify  the  subscene  he  is  interested  in  and  to 
write  it  onto  another  tape.  This  increases  tape  processing 
efficiency  and  reduces  computation  cost. 

SUBGM  Version  of  SUBSET  that  is  used  to  geometrically  correct,  rotate,  and 
scale  Landsat  data  during  extraction  of  a  subject. 

NMAP  Creates  a  printed  grey  scale  brightness  map  which  shows  the  overall 
scene  (or  SUBSET)  features.  The  program  assists  the  user  in 
recognizing  and  visually  correlating  the  data  on  tape  to,  either 
SUBSET  or  SUBGM,  with  areas  seen  on  photographic  imagery  or  on 
topographic  maps. 

STATS  Computes  and  outputs  statistics  such  as  the  means  and  standard 

deviation  of  the  digital  numbers,  correlation  matrices,  and  frequency 
histograms  for  specified  areas  and  channels. 

STATS2  has  in  addition  the  capability  of  producing  two  dimensional 
histograms  of  the  responses  of  two  spectral  channels  has  been  added. 
PPD  Classifies  data  using  parallelpiped  algorithm.  The  spectral 

characteristics  of  classes  are  defined  by  upper  and  lower  bounds  for 
each  channel. 

CLUS  Unsupervised  classifier  which  develops  a  set  of  spectral  signatures 
by  clustering. 

CLASS  Performs  supervised  classification  using  the  Euclidean  distance 

algorithm.  The  CLASS  routine  was  modified  to  allow  the  output  of  raw 
Landsat  data  (digital  numbers)  to  disk  or  tape.  This  version  is 
called  CLDISK. 

Both  the  STATS2  and  the  CLDISK  programs  are  only  available  on  the  USDA 
computer.  Correction,  calibration  and  further  analysis  of  the  Landsat  data  is 
performed  on  the  WCC’s  IBM  370  computer  using  the  Statistical  Analysis  System 
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(SAS)  of  SAS  Institute  Inc.  (Helwig  and  Council,  1979).  Examples  of  the  above 
programs  and  their  outputs  are  given  in  the  Appendices  2  and  3* 
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3.  Preparation  of  data 

3.1  Curve  numbers  and  watersheds 

Rainfall  and  runoff  records  were  obtained  for  90  watersheds  in  the 
western  USA  from  various  sources,  such  as  the  United  States  Geological  Survey 
(USGS) ,  United  States  Forest  Service  (USFS),  USDA-ARS,  and  others.  The  USDA- 
ARS  data  were  retrieved  from  the  Water  Data  Bank  on  the  USDA  Washington 
Computer  Center  (WCC)  utilizing  the  REPHLEX  procedures  described  by  Thurman, 
et  al.  (1983)*  For  each  watershed  the  curve  numbers  were  computed  using  three 
different  methods  and  two  input  data  arrangements.  Thus  a  total  of  six  CN 
were  determined  per  watershed.  The  computational  procedures  for  each  method 
are  given  below. 

Mean  Curve  numbers  for  all  the  events  are  computed  and  the 


arithmetic  mean  is  taken. 


Non-Linear  Equation  5  and  6  of  section  2.2  result  in  the 


regression  relationship. 


Q  =  (P+2-200/CN)2/(P-8+8.00/CN) ,  P  >  (200/CN)-2 


[133 


Q  =  0 


,  P  <  ( 200/CN) -2 


Equation  13  is  used  in  an  algorithm  that  fits  the  general 


model  through  all  the  points  and  finds  the  optimum  CN,  i.e.  the 


CN  for  which  the  sum  of  squares  of  deviations  is  smallest,  or 

p 

where  (Q-Q)  has  a  minimum. 


Relative 


Found  by  trial-and-error  on  the  computer,  using  the 


storm  size  following  procedure. 


1 .  Calculate  average  S  of  the  events  and  take  the  smallest 


rainstorm. 
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2.  Check  if  it  meets  the  criterion  P/S  >  0.46. 

3  V 

3.  If  P/S„,.  <  0.46,  drop  this  event,  and  go  back  to 

3  V  — 

step  1. 

4.  Take  the  mean  CN  of  the  events  left  in  the  data  set. 

These  computations  were  carried  out  for  90  watersheds,  and  the  results 

are  presented  in  Appendix  4.  Sample  SAS  and  FORTRAN  programs  to  prepare  and 
analyse  the  hydrologic  data  are  displayed  in  Appendix  5. 

Unfortunately,  for  1 6%  of  the  watersheds  the  selection  method  did  not 
come  up  with  a  CN  for  the  ordered  P  and  Q  arrays.  For  the  remaining  70 
watersheds  which  yielded  a  CN,  the  CN  using  the  "mean”  method  is  on  the 
average  3-7  CN  higher  than  the  "selection"  method,  for  the  ordered  data  sets, 
while  the  "regression"  method  was  only  0.2  CN  higher.  Thus  although  the  CN 
computation  using  only  large  storms  is  theoretically  more  correct  than  a  least 
squares  fit,  in  practice  there  is  very  litle  difference  between  the  results  of 
the  two  methods.  This  is  why  the  non-linear  regression  method  is  selected  for 
further  use  in  this  study.  Table  4  shows  the  average  CN  for  the  different 
methods,  including  only  the  70  watersheds  for  which  a  CN  could  be  calculated 
using  the  selection  method  on  natural  and  ordered  events.  Figure  10 
visualizes  the  differences  between  the  various  average  CN,  and  Figure  11 
illustrates  the  CN  found  by  non-linear  regression  for  Chickasha  R5  using  both 
the  natural  and  the  ordered  events. 


Table  4.  Average  curve  numbers  for  different  calculation  methods 


Method 

Data  arrangement 

Natural 

Ordered 

Mean 

83.8 

83.0 

Non-linear  regression 

78.6 

79.5 

Relative  storm  size 

77.9 

79.3 
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A  selection  of  basins  for  the  initial  study  was  made  from  all  90 
watersheds  available  based  on  the  following  criteria  : 

1.  the  size  has  to  be  at  least  20  acres  (preferably  50  acres)  to  facilitate 
location  on  Landsat  tapes. 

2.  there  should  be  a  wide  range  in  CN. 

3.  variation  in  land  cover  type. 

4.  location  in  various  parts  of  the  West. 

5.  Landsat  tapes  should  be  readily  available  for  immediate  analysis  (i.e.,  no 
long  waits  for  tape  orders). 


NATURAL 


ORDERED 


82.8 


Mean 


83.0 


Figure  10.  Results  of  the  different  CN  methods. 
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Figure  11.  Non  Linear  regression  for  Chickasha  R5  data  using  both  natural  and 
ordered  P-Q  data. 

The  watersheds  meeting  these  criteria  are  listed  in  Table  5.  The 
location  of  the  watersheds  is  pictured  in  Figure  12. 

3.2  Landsat  response 

3.2.1  Sources 

Landsat  images  were  obtained  from  the  tape  library  of  the  Remote  Sensing 
Systems  Laboratory  (RSSL)  at  the  University  of  Maryland's  Department  of  Civil 
Engineering  and  from  files  at  the  Hydrology  Laboratory.  The  scenes  were  all 
originally  purchased  from  the  Earth  Resources  Observation  Systems  (EROS)  data 
center  in  Sioux  Falls,  SD.  Table  6  identifies  the  tapes  used  in  the  analysis. 


Table  5.  Summary  of  Watershed  Characteristics 


CO 

p 

cO 

c 

43- 

X- 

X- 

X- 

X- 

p 

p 

CM 

m 

on 

•o 

o 

X- 

c — 

e'¬ 

X- 

X- 

p 

p 

X- 

co 

X- 

•rH 

cx 

cx 

er. 

cr 

CX 

cx 

cx 

cr 

CX 

cx 

Q— f 

4-> 

r — 

T — 

, — 

, — 

T - 

, — 

, — 

t — 

• — 

t — 

O 

o 

1 

1 

i 

l 

1 

1 

1 

1 

1 

1 

CP 

CX 

X- 

x- 

X- 

x- 

0 

O 

P 

CM 

p 

•D 

1 — 1 

p 

P 

P 

p 

p 

43- 

43- 

P 

CO 

p 

O 

1 — 1 

cx 

C x 

cx 

cx 

CX 

cx 

cx 

CX 

CX 

cx, 

•H 

o 

, — 

, — 

t — 

, — 

, — 

, — 

T - 

t — 

t — 

t — 

P 

o 

CP 

(X 

Cl-, 

o 

CO 

CO 

co 

CO 

CO 

co 

co 

co 

CO 

co 

(X 

(X 

CX 

cx 

cx 

Cx-I 

x 

a 

0 

0 

CP 

CO 

«=c 

<c 

< 

eC 

*% 

co 

CO 

co 

co 

co 

a 

p 

DP 

X 

X 

X 

X 

p 

CO 

3 

•o 

o 

CO 

2 

cx 

43- 

43- 

P 

0 

■xr 

CM 

on 

X- 

0 

O 

X- 

X- 

OO 

CO 

P 

p 

cx 

00 

cx 

CO 

CO 

• 

• 

CO 

•o 

■D 

CO 

P 

c 

C 

CO 

bO 

o 

O 

co 

co 

o 

O 

p 

CO 

• 

• 

bO 

CO 

co 

■o 

TJ 

•D 

T3 

p 

p 

o 

O 

C 

C 

•• 

CP 

bC 

o 

O 

O 

O 

P 

Cm 

bO 

bC 

O 

O 

CO 

•rH 

•> 

3 

c 

CO 

•> 

•• 

P 

P 

P 

0 

co 

CO 

CO 

(X 

P 

CO 

CO 

O 

•H 

P 

0 

P 

p 

p 

w 

3 

CO 

CO 

O 

CO 

3 

3 

3 

> 

p 

co 

CO 

a 

Cm 

•> 

•> 

P 

P 

p 

o 

P 

P 

P 

CO 

CO 

P 

P 

p 

o 

CO 

bD 

bD 

— 

•> 

T3 

TJ 

CO 

CO 

CO 

CP 

CP 

O 

O 

p 

CP 

CP 

p 

P 

O 

O 

p 

p 

p 

P 

> 

> 

3 

3 

2 

5 

p 

p 

p 

CP 

•rH 

•rH 

P 

P 

TJ 

■O 

CP 

CP 

a) 

CO 

pi 

P 

CO 

CO 

p 

P 

CO 

CO 

(O 

CP 

CO 

CO 

CO 

cO 

CO 

CO 

CP 

CP 

<p 

Q 

2 

z 

P-I 

cx 

X 

X 

0 

Q 

Q 

o 

x- 

CM 

CM 

p 

CO 

• 

• 

• 

• 

_ _ 

o 

x- 

X- 

cx 

X- 

43- 

p 

CO 

p 

CM 

CP 

P 

C\J 

CM 

1 — 

CM 

co 

p 

T - 

0 

43- 

N 

P 

43- 

, — 

p 

p 

43- 

•rH 

co 

CP 

p 

cO 

CO 

bd 

CX 

IX 

bC 

H 

Sh 

>H 

Sh 

>H 

P 

<c 

o 

O 

0 

O 

b 

3 

CO 

<P 

C 

c 

CP 

>> 

43- 

in 

P 

X- 

CO 

a> 

<p 

2 

cx 

CX 

cx 

CX 

CP 

CP 

CP 

p 

p 

CP 

(P 

2 

0 

a 

CP 

co 

CO 

CO 

CO 

0 

p 

p 

CO 

w 

c 

P 

p 

P 

X 

0 

0 

P 

>> 

DC 

o 

CO 

co 

CO 

CO 

>> 

O 

p 

CO 

p> 

CO 

co 

CO 

CO 

CO 

CO 

p 

a 

cc 

CO 

X 

5 

•H 

3 

w 

P 

o 

0 

V 

0 

C*-H 

P 

O 

c 

• 

E-h 

E 

•H 

•H 

•H 

•rH 

rH 

P 

bO 

CO 

ClH 

< 

o 

p 

P 

P 

X 

CO 

O 

3 

p 

• 

2 

E-* 

o 

O 

O 

0 

X 

S 

Q 

pH 

:s 

-z  o- 


Figure  12.  Location  of  the  watersheds  selected. 

3.2.2  Mapping  and  data  retrieval  procedures 

The  watersheds  that  have  been  selected  in  section  3»1  were  located  on 
USGS  1:24,000  and  1:250,000  scale  topographic  maps  and  (if  available)  on 
aerial  photographs.  For  an  initial  look  at  the  scenes  use  was  made  of  an 
image  display  system  at  the  RSSL.  The  image  processing  system  utilized  was  an 
International  Imaging  System,  model  70F,  interfaced  to  a  Digital  PDP  11/24  RSX 
-  1 1 M  operating  system.  The  results  of  processing  were  displayed  on  a 
Mitsubishi  512  x  512  pixel  color  monitor.  Comparing  the  color  display  and  the 
topographic  maps,  the  watershed  boundaries  were  roughly  described  in  terms  of 
their  line  and  element  numbers.  Using  this  information  from  the  RSSL  system  a 
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Table  6.  Tape  information 


Watershed 

Path-Row 

Scene  ID 

Date  of 

observation 

Tombstone  W4 

38-38 

10984-17093 

03  APR  75 

Halfway  and  Morris 

41-32 

20155-17310 

26  JUN  75 

Dugout  Creek 

38-30 

11416-16394 

08  JUN  76 

Frank  Draw 

37-30 

11073-16554 

01  JUL  75 

West  Fork  Dry  Cheyenne 

39-30 

20153-17184 

24  JUN  75 

Chcikasha  R5,  R6,  R7  and  R8 

30-36 

1256-16415 

05  APR  73 

1400-16402 

27  AUG  73 

1508-16383 

13  DEC  73 

1706-16332 

29  JUN  74 

1814-16295 

15  OCT  74 

SUBSET  was  made  of  the  general  area  and  was  visualized  by  a  line  printer  map 
(NMAP-  procedure).  Subsequently  a  PPD  program,  using  the  rather  unique 
response  pattern  of  open  water,  was  performed  on  the  same  SUBSET  area,  thus 
creating  a  line  printer  map  that  only  indicates  lakes,  farm  ponds,  etc.  This 
map  plus  the  detailed  information  from  the  topographic  maps  and  aerial 
photographs  was  sufficient  to  accurately  indicate  watershed  boundaries  on  the 
line  printer  map.  An  example  of  the  location  of  a  watershed  is  presented  in 
Appendix  6.  The  digital  numbers  (DN)  within  the  watershed  boundaries  were 
written  to  disk  files  by  means  of  the  CLDISK  program. 
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3.2.3  Correcting  and  calibrating  procedures. 

Landsat  digital  images  are  commonly  analyzed  by  using  the  digital  numbers 
(DN)  for  each  pixel.  Although  this  procedure  may  be  satisfactory  when  only  a 
single,  internally  consistent  image  is  used,  in  the  case  of  this  study  the 
procedure  would  produce  incorrect  results,  since  the  DN  were  collected  by 
different  satellites  and  at  different  times.  Therefore,  the  DN  should  be 
converted  to  dimensioned  equivalents  such  as  radiance  and  reflectance.  The 
digital  levels  are  related  by  a  linear  model  (calibration)  to  the  intensity  of 
reflected  radiant  energy  that  are  not  directly  comparable  among  the  Landsat 
satellites  because  of  differences  in  the  calibration  of  their  multispectral 
scanner  (MSS)  instruments.  In  addition,  comparison  or  combinations  of  images 
taken  at  different  times  requires  correction  for  different  angles  of  solar 
illumination. 

The  radiance  as  measured  at  the  satellite  in  a  single  band  is  calculated 
by 

Radiance  =  DN/Dmax  (Lmax  -  Lmin)  +  Lmin  (mW.cm2“sr"1 )  [143 

where 

DN  =  the  digital  value  of  a  pixel  from  the  CCT  (-) 

Dmax  =  the  maximum  digital  value  of  that  band  (-) 

Lmax  =  radiance  measured  at  detector  saturation  (mWcm-^  sr-1) 

Lmin  =  the  lowest  radiance  measured  by  detector  (mWcm-2  sr”1) 

Table  7  gives  the  Lmax,  Lmin,  and  Dmax  values  for  Landsats  1,  2,  and  3, 
and  for  different  scanner  calibrations  at  different  times.  It  is  modified 
from  the  Landsat  Data  Users  Handbook  (USGS,  1979,  p  AE-16)  and  Robinove  et  al. 
(1981).  Although  no  Landsat-3  data  were  included  in  this  analysis,  the  values 


are  included  for  reference. 


Table  7.  Landsats  1,  2,  and  3  detector  response. 
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The  reflectance,  which  is  the  percentage  of  radiance  to  irradiance,  is 
calculated  in  a  single  band  for  a  Lambertian  surface  by 


Reflectance  =  71 /E  sin  a 


DN/Dmax(Lmax-Lmin)+Lmin 


(  -  ) 


[15] 


where 


a  =  solar  elevation,  measured  from  the  horizontal,  as  annotated  on 
Landsat  images. 

E  =  average  solar  irradiance  in  mW  cm-^  at  the  top  of  the 
atmosphere : 

Band  4  =  17*70  mW  cm-2 
Band  5  =  15.15  mW  cm-2 
Band  6  =  12.37  mW  cm-2 
Band  7  =  24.91  mW  cm-2 

Solar  irradiance  may  vary  as  much  as  7%  from  aphelion  to  perihelion.  However, 
for  the  normal  range  of  surface  reflectance,  the  variation  in  reflectance 
would  be  less  than  1  % .  Other  corrections,  for  atmospheric  absorption, 
atmospheric  scattering,  and  the  radiance  effects  of  nearby  pixels  (adjacency 
effects)  have  not  been  applied.  Corrections  for  these  factors  could  be  done 
in  a  highly  sophisticated  manner,  such  as  the  use  of  radiosonde  observations 
at  the  time  of  imaging,  in  order  to  determine  the  state  of  the  atmosphere. 

Because  these  required  measurments  were  not  available,  atmospheric 
effects  were  not  taken  into  account. 

Four  assumptions  were  made  in  the  reflectance  calculations: 

1.  the  state  of  the  atmosphere  is  the  same  for  all  scenes 

2.  the  terrain  surface  is  a  Lambertian  reflector 

3.  the  average  terrain  slope  is  zero 

4.  the  sun  angle  contribution  to  the  scene  brightness  is  uniform  over  the 


entire  scene 
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The  values  of  Table  7  in  equation  T.153  give  the  following  equations  used 
to  compute  the  reflectance  for  each  pixel: 

Landsat  1 


Band  4 
Band  5 
Band  6 
Band  7 


Reflectance  =  (DN/289)/  sin  a 
=  (DN/306)/  sin  a 
=  (DN/284)/  sin  & 
=  (DN/125) /  sin  a 


20 

21 

22 


23j 


Landsat  2,  before  7/16/75 

Band  4  Reflectance  =  0.1775  ( (DN/63.50)+0. 10)/sin  a 

Band  5  =  0.2074  ( (DN/85 .81 )+0 .07) /sin  a 

Band  6  =  0.2540  ((DN/95.49)+0.07)/sin  a 

Band  7  =  0.1261  ( (DN/15.71 )+0.l4)/sin  a 

From  the  above  calculations  spectral  reflectance  patterns  were 

constructed  for  the  ten  watersheds  that  were  selected  in  section  3»1  using  the 
June  74  scene  for  the  Chickasha  watersheds.  These  patterns  are  shown  in 
Figure  13» 

To  illustrate  that  there  is  a  probability  associated  with  the  response  in 
each  band,  Figure  14  gives  the  distribution  histograms  in  band  5  for  three 
watersheds . 


3.2.4  Reflectance  index  models 

Variation  in  CNs  between  the  various  watersheds  is  caused  by  differences 
in  factors  such  as  1)  hydrologic  condition  of  the  surface;  2)  soil 
characteristics;  3)  climate;  and  4)  topography.  These  factors  can  not  be 
measured  directly  by  a  remote  sensing  system.  However,  remote  sensing  tends 
to  give  an  integrated  result  of  the  factors  mentioned  above. 
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Figure  13.  Spectral  reflectance  patterns  for  10  watersheds. 


The  link  between  Landsat  reflectance  and  CNs  is  made  by  presenting  the 
multiband  data  in  a  single  value,  the  reflectance  index.  The  reflectance 
index  is  a  quantity  obtained  directly  or  by  ratioing,  differencing,  or 
otherwise  transforming  spectral  data  in  a  Reflectance  Index  Model  (RIM),  to 
represent  plant  canopy  characteristics  such  as  leaf  area  index,  biomass,  and 
percent  cover,  and  soil  properties.  In  the  literature  numerous  RIMs  are 
described  and  tested  (Richardson  and  Wiegand,  1977;  Jackson  et  al.,  1980; 
Perry  and  Lautenschlager ,  1984),  and  in  this  section  a  limited  number  of  them 
will  be  summarized  and  employed  in  the  CN  study. 
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Figure  14.  Comparison  of  the  reflectance  distributions  in  band  5  for  three 
watersheds 


The  actual  reflectance  values  from  the  KSS  channels  are  used  in  the 
computations  rather  than  the  digital  counts  which  can  not  be  compared  between 
the  different  scenes.  Most  formulae  are  based  on  ratios  or  linear 
combinations  and  exploit  differences  in  the  reflectance  patterns  of  green 
vegetation  and  other  objects  as  summarized  in  Figure  15. 

Reflectance  values  from  the  individual  channels  (CH4,  CH5,  CH6,  and  CH7) 
have  been  used  previously  to  estimate  percent  ground  cover  and  vegetative 
biomass. 

The  ratio  of  reflectance  values  from  two  bands  is  a  simple  and  useful 
vegetation  index,  if  the  bands  are  properly  chosen.  One  criterion  for 
choosing  two  bands  for  a  ratio  vegetation  index  is  that  the  data  from  one  band 
should  decrease  with  increasing  green  vegetation,  and  data  from  the  other  band 
should  increase  with  increasing  green  vegetation.  This  is  the  case  for  the 
red  band  and  the  near  infrared  bands,  as  can  be  seen  from  figure  15. 
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Figure  15.  Idealized  reflectance  patterns  of  vegetation  and  soil. 


In  this  study  the  following  ratios  are  applied: 


R54  =  CH5/CH4 

[24] 

R65  *=  Cl}6/CH5 

[253 

R75  =  CH7/CH5 

[263 

R65  and  R75  have  been  shown  to  be  sensitive  to  the  amount  of  vegetation 
whereas  R54  distinguishes  best  between  soil  and  rock  groups  (Tucker,  1979; 
Eliason  et  al.,  1981). 

The  ratio  of  the  difference  between  values  for  two  bands  and  the  sum  of 
the  values  for  the  two  bands,  the  normalized  difference  (ND),  is  used  for 
monitoring  vegetation.  This  model  was  developed  by  Deering  and  Haas  (1973) 
and  discussed  by  Rouse  et  al.  (1973).  Deering  et  al.  (1975)  added  0.5  to  ND 
to  avoid  negative  values  and  took  the  square  root  of  the  result  to  stabilize 
the  variance.  This  index  is  referred  to  as  the  transformed  vegetation  index 
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and  will  be  denoted  by  TVI.  Both  ND  and  TV I  were  proposed  using  CH7 ,  CH6 ,  and 
CH5: 


ND6  =  ( CH6  -  CH5)/(CH6  + 

CH5) 

L27  3 

ND7  =  (CH7  -  CH5)/(CH7  + 

CH5) 

[28] 

TVI6  =  (ND6  +  0.5)1/2 

[293 

TVI7  =  (ND7  +  0.5) 1/2 

[303 

In  addition  to  ratios  and  difference  ratios,  many  other  combinations  of 
spectral  bands  have  been  used  as  reflectance  indices.  Kauth  and  Thomas  (1976) 
have  developed  a  technique,  using  vector  analysis  in  four  dimensional  space  to 
produce  an  orthogonal  transformation  of  the  original  Landsat  data  space.  They 
named  the  four  new  axes  Soil  Brightness  Index  (SBI),  Green  Vegetative  Index 
(GVI),  Yellow  Vegetative  Index  (YVI)  and  non  such  (NSI): 


SBI  =  0.332  CH4  +  0.603CH5  +  0.675CH6  +  0.262  CH7  H 3 1  □ 
GVI  =  -0.283  CH4  -  0.660CH5  +  0.577CH6  +  0.388  CH7  C 32l 
YVI  =  -0.899  CH4  +  0.428CH5  +  0.076CH6  -  0.041  CH7  C333 
NSI  r  -0.016  CH4  +  0.131CH5  -  0.452CH6  +  0.882  CH7  C 34D 


Richardson  and  Wiegand  (1977)  used  the  perpendicular  distance  to  the  "soil 
line"  as  an  indicator  of  plant  development.  The  "soil  line,"  a  two 
dimensional  analogue  of  the  Kauth-Thomas  SBI,  was  estimated  by  linear 
regression.  Two  perpendicular  vegetation  indices  were  proposed: 

PVI6  =  [ (-2.507-0.457CH5  +  0.498CH6)2 
+  (2.734  +  0.498CH5  -  0.543CH6)2J 1/2  II35] 

PVI7  =  [(0.355CH7  -  0.149CH5)2 

+  (0.355CH5  -  0.852CH7)2:172  [36] 

The  difference  Vegetation  Index  (DVI),  suggested  as  computationally 
easier  than  PVI7,  is  essentially  a  rescaling  of  PVI7  (Richardson  and  Wiegand, 
1977): 
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DVI  =  2.4  CH7  -  CH5  H 371 

The  Ashburn  Vegetation  Index  (AVI)  was  suggested  by  Asburn  (1978)  as  a 
measure  of  green  growing  vegetation: 

AVI  =  2.0CH7  -  CH5  l  381 

3.2.5  Reflectance  index  models  and  ground  truth. 

Ground  truth  was  collected  at  approximately  the  same  time  as  the  Landsat 
overpass  for  four  of  the  five  available  scenes  of  the  rangeland  watersheds  in 
Chickasha,  Oklahoma.  In  the  watersheds  (R5,  R6,  R7,  and  R8)  the  following 
information  was  obtained:  1)  a  field  estimate  of  the  standing  vegetation  that 
is  green  matter;  2)  fresh  biomass  (in  grams  per  sample  plot);  and  3)  dry 
biomass  (in  gm  per  sample  plot). 

Correlations  (r)  were  made  of  ground  truth,  acquired  from  four  watersheds  at 
four  different  times  (N  =  16),  with  the  19  RIMs  which  are  discussed  above. 

The  results  are  arranged  in  Table  8. 

In  general  the  percentage  of  green  matter  shows  the  highest  correlations. 
Those  RIMs  using  band  5  and  6  combinations  correlate  better  with  the  ground 
truth  information  than  the  models  that  use  band  6  and  7  combinations. 

Similar  results  were  found  by  Richardson  and  Wiegand  (1977). 

Perry  and  Lautenschlager  (1984)  demonstrated  that  R65,  ND6  and  TVI6  are 
equivalent  for  decision  making  (and  so  are  R75,  ND7  and  TVI7).  This  finding 
is  again  underlined  by  the  closeness  of  the  correlation  coefficient  for  these 
RIMs. 

Surprising  is  the  high  negative  correlation  of  PVI6. 
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Table  8. 

Correlations  between 

reflectance  index  models  and 

ground  truth 

RIM 

Green  matter  (%) 

Fresh  biomass 

dry  biomass 

CH4 

0.211 

0.186 

0.153 

CH5 

-0.310 

-0.416 

-0.441 

CH6 

0.375 

-0.002 

-0.046 

CH7 

0.369 

0.001 

-0.003 

R54 

-0.480 

-0.581** 

-0.591** 

R65 

0.780* 

0.557** 

0.551** 

R75 

0.619** 

0.449 

0.477 

ND6 

0.796* 

0.545** 

0.534** 

ND7 

0.626* 

0.461 

0.483 

TVI6 

0.801* 

0.541** 

0.529** 

TVI7 

0.627* 

0.463 

0.484 

GVI 

0.803* 

0.396 

0.388 

SB  I 

0.117 

-0.157 

-0.190 

YVI 

-0.518 

-0.648* 

-0.652* 

NS  I 

-0.442 

-0.050 

-0.034 

PVI6 

-0.801* 

-0.293 

-0.254 

PVI7 

0.594** 

0.226 

0.234 

DVI 

0.594** 

0.226 

0.231* 

AVI 

0.627* 

0.271 

0.281 

*  Statistically  significant  at  0.01  probability  level 

**  Statistically  significant  at  0.05  probability  level 
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Results 

4.1  Relationships  between  curve  numbers  and  reflectance  index  models. 

First  a  linear  relationship  was  tested  for  the  curve  numbers  and 
reflectance  indices,  which  were  derived  as  is  explained  in  the  previous 
chapter.  Figure  16  shows  the  plots  for  R65  and  SBI. 


CN 


i  0  <2  1  *  <6  16  20  22 

r©s 


CN 


Figure  16.  Linear  relationships  for  curve  numbers  and  two  RIMs  (N=10). 

Correlations  for  all  19  RIMs  and  CNs  were  established  and  are  given  in 
Table  9. 

Notable  are  the  high  values  for  the  single  channels  4,  5,  and  6.  Again 
the  RIMs  that  consist  of  band  6  and  5  combinations  give  higher  correlations 
than  the  band  7  and  5  combinatons.  If  we  take  only  the  band  6  and  5 
combinations  and  the  RIMs  that  show  a  correlation  which  is  statistically 
significant  at  the  0.01  probability  level,  five  RIMs  remain.  Descriptions  of 
the  relationships  between  these  RIMs  and  CNs  were  established  using  linear 
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regression.  The  parameters  (al  and  bl)  of  the  general  relationship 

CN  =  a 1  +  bl  (RIM) 

and  the  correlation  coefficient  (r)  belonging  to  it  are  given  in  Table 
10.  Also  shown  in  table  10  are  the  parameters  and  correlations  for  the 
general  equation 

In  CN  =  a2  +  b2  (RIM) 

which  suggests  the  exponential  relationship 
CN  =  a3  x  eb2 


Table  9.  Correlations  between  curve  numbers  and  reflectance  index  models 


(N  =  10) 


RIM 

r 

RIM 

r 

CH4 

0.854* 

TVI6 

-0.841* 

CH5 

0.899** 

TVI7 

-0.808* 

CH6 

0.851* 

SBI 

0.918** 

CH7 

-0.095 

GVI 

-0.705 

R54 

0.636 

YVI 

-0.674 

R65 

-0.865* 

NS  I 

-0.077 

R75 

-0.840* 

PVI6 

0.685 

ND6 

-0.850* 

PVI7 

-0.524 

ND7 

-0.821* 

DVI 

-0.524 

AVI 

-0.571 

*  Statistically  significant 

at 

0.01  probability  level 

**  Statistically 

significant 

at 

0.001  probability  level 

Table  10  indicates  that  an  exponential  relationship  generally 

improves  the 

fit,  but  probably  not  significantly. 


Personal  communication  with  Ken  Renard,  director  of  the  Southwest 
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Rangeland  Watershed  Research  Center  in  Tucson,  (Arizona)  cleared  up  the  poor 
fit  of  the  Tombstone  watershed  (CN  =  79  in  Figure  16).  This  watershed  has  a 
very  high  infiltration  in  the  main  channel.  Thus  it  has  generally  a  low 
discharge  at  the  flow  measurement  structure  and  therefore  an  underestimated 


Table  10.  Parameters  of  linear  and  exponential  relationships  between  curve 
numbers  and  reflectance  index  models. 


RIM 

linear 

exponential 

al 

bl 

r 

a2 

b2 

r 

CH5 

47.55 

212.19 

0.899 

3.940 

2.840 

0.915 

R65 

109.56 

-20.72 

-0.865 

4.774 

-0.280 

-0.888 

ND6 

91.33 

-71.27 

-0.850 

4.528 

-0.960 

-0.870 

TVI6 

176.15 

-118.85 

-0.841 

5.669 

-1.599 

-0.860 

SB  I 

17.41 

183.20 

0.918 

3.537 

2.453 

0.934 

CN .  The 

actual  CN  on 

the  slopes  would  be  an  estimated 

10-15  CN  higher. 

Leaving 

this  watershed 

out  the 

data  set,  gives 

the  following  correlation 

coefficients . 

Table  11 

.  Correlations  between 

curve  numbers 

and  reflectance  index 

models 

(N  =  9) 

RIM 

linear 

exponential 

CH5 

0.967 

0.976 

R65 

-0.927 

-0.944 

ND6 

-0.950 

-0.962 

TVI6 

-0.955 

-0.965 

SBI 

0.938 

0.950 
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The  parameters  of  both  the  linear  and  exponential  relationships  are  presented 
in  Appendix  7. 

4.2  Significance  of  seasonal  fluctuation  of  Landsat  data. 

Figure  17  illustrates  the  seasonal  fluctuation  of  the  R65  and  SBI  for  the 
Chickasha  watershed  R7.  To  find  out  what  the  effect  of  this  seasonality  is  on 
the  relationships  shown  in  Table  10,  R65  and  SBI  values  from  the  Chickasha 


R6S 


DATE 
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DATE 

CHICKASHA  R7 


Figure  17.  Seasonal  fluctuation  of  R65  and  SBI. 

watersheds  were  calculated  for  five  different  times  of  the  year  and 
successively  used  in  the  linear  regression.  The  correlation  coefficients  were 
calculated  and  are  presented  in  Table  12. 
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Table  12. 

Correlation 

coefficients 

for  two 

reflectance  index  models 

and  curve 

numbers  (N  = 

10) . 

RIM 

April 

* 

July 

August 

October 

December 

R65 

-0.853 

-0.865 

-0.880 

-0.896 

-0.793 

SB  I 

0.881 

0.918 

0.913 

0.841 

0.898 

*  from  Table  10. 

-  Data  from  Chickasha  watersheds  vary  over  the  year. 

-  Data  from  other  watersheds  are  from  April  through  June. 

In  Table  12  it  can  be  seen  that  the  use  of  reflectance  data  (in  grassland 
watersheds)  from  different  times  of  year  for  the  Chickasha  watersheds  do  not 
significantly  change  the  correlation  of  a  linear  relationship.  This  indicates 
that  the  seasonal  fluctuation  of  these  RIMs  could  be  relatively  small  as 
compared  to  the  differences  between  the  actual  rangeland  classes. 
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5.  Conclusions  and  recommendations. 

1 .  High  correlations  between  rangeland  runoff  curve  numbers  and  various 
reflectance  indices  obtained  from  Landsat  MSS  data  resulted  through 
regression  analysis.  Logarithmic  transformations  slightly  improved  the 
relationships.  Discarding  an  anomalous  basin  in  Arizona  significantly 
improved  the  regression  relationships. 

2.  Five  reflectance  indices  were  selected  from  a  total  of  19  originally 
calculated  from  the  Landsat  data  because  of  their  high  correlations.  Of 
these  five  indices,  R65,  ND6 ,  and  TVI6  provide  virtually  identical 
information  on  vegetation  conditions.  The  other  two  indices,  CH5  and  SBI , 
provide  information  on  soil  background. 

3.  Although  the  "selection"  method  for  determining  CNs  was  preferred  over  the 
mean  and  regression  methods,  the  regression  method  was  chosen  because  it 
provided  nearly  similar  results  to  the  "selection"  method  and  it  is  a  lot 
less  effort.  Ordered  data  are  used  rather  than  the  natural  data  because 
they  result  in  a  more  stable  CN  value. 

4.  Future  work  should  test  the  applicability  of  this  method  over  a  wide*” 
variety  of  rangeland  conditions. 

-Test  more  basins. 

-Determine  best  season,  if  there  is  one. 

-Attempt  to  see  if  a  combination  of  RIMs  would  provide  even  better  CN 
estimation . 

-Evaluate  other  uses  of  this  approach  to  SCS  needs. 

-Based  on  all  this,  come  up  with  a  final  methodology,  of  which  an  example 
is  given  in  Figure  18. 


Figure  18.  Outline  for  a  computerized  curve  number  estimation  for  rangeland. 

A  Geographic  Information  System  (GIS)  has  the  capability  of  integrating 
different  sources  of  information  related  to  geographicaly  the  same 
area.  Various  GIS’s  are  available  on  the  market.  The  spectral  analysis 
program  should  incorporate  the  CN-RIM  relationships. 
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List  of 

AMC 

AVI 

CCT 

CH4 

CH5 

CH6 

CH7 

CN 

DN 

Dmax 

E 

F 

GVI 

la 

Lmax 

Lmin 

MSS 

N 

NEH-4 

ND6 

ND7 

NS  I 

ORSER 

P 

PVI6 

PVI7 

Q 
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RIM 

R54 

R65 

R75 

S 

SAS 
SB  I 
SCS 
TVI6 
TV  17 
USDA 
USGS 
WCC 
YVI 
O'. 


abbreviations  and  symbols 


antecedent  moisture  class 
Ashburn  vegetation  index 
computer  compatible  tape 
reflectance  in  channel  4 
reflectance  in  channel  5 
reflectance  in  channel  6 
reflectance  in  channel  7 
runoff  curve  number 
digital  number 

maximum  digital  value  in  a  band 

average  Solar  irradiance  at  top  of  atmosphere  (mWcm-^) 
total  acutal  retention  (in.) 
green  vegetative  index 
initial  abstraction  (in.) 

radiance  measured  at  detector  saturation  (mWcnT^Sr  ) 
lowest  radiance  measured  by  detector  (mWcm  Sr  ) 
multispectral  scanner 
number  of  data  points 

national  engineering  handbook,  section  4 

normalized  difference,  using  bands  5  and  6 

normalized  difference,  using  bands  5  and  7 

"non  such"  vegetation  index 

office  for  remote  sensing  of  earth  resources 

precipitation  (in.) 

perpendicular  vegetation  index 

perpendicular  vegetation  index 

surface  runoff  (in.) 

correlation  coefficient 

reflectance  index  model 

ratio  of  channel  5  over  channel  4 

ratio  of  channel  6  over  channel  5 

ratio  of  channel  7  over  channel  5 

maximum  retention  (in.) 

statistical  analysis  system 

soil  brightness  index 

Soil  Conservation  Service 

transformed  vegetation  index 

transformed  vegetation  index 

United  States  Department  of  Agriculture 

United  States  Geological  Survey 

Washington  Computer  Center 

yellow  vegetation  index 

sun  elevation  angle 


-49- 


Appendices 

1.  Runoff  probability,  relative  storm  depth,  and  curve  numbers. 

2.  Examples  of  Job  Control  Language  for  the  ORSER  programs  on  the  IBM  370 
computer,  and  their  outputs. 

3.  Examples  of  SAS  programs  used  in  the  Landsat  analysis. 

4.  Summary  of  curve  number  computations. 

5.  Examples  of  some  SAS  and  FORTRAN  programs  used  in  the  curve  number 
analysis. 

6.  Location  of  Tomstone,  W4  on  line  printer  map. 

7.  Parameters  of  linear  and  exponential  relationships  between  curve  numbers 
and  reflectance  index  models  for  9  watersheds. 


AP-1 


Appendix  1.  Runoff  Probability,  Relative  Storm  Depth,  and  Curve  Numbers 
Richard  H.  Hawkins^-,  M.  ASCE,  Allen  T.  Hjelmfelt,  Jr. 2  M.  ASCE, 

and  Adrian  W.  Zevenbergen2 

Using  existing  antecedent  moisture  condition  ( AMC) — Curve  Number 
relationships,  a  general  expression  of  rainf al 1 -runoff  was  developed  for  AMC  I 
and  AMC  III.  The  distribution  of  the  probability  of  event  runoff  exceeding 
zero  is  found  to  be  approximately  lognormal.  Relative  storm  size  is  then 
proposed  to  be  defined  on  the  ratio  P/S,  where  a  "large"  storm  has  P/S>0.46, 
when  90  percent  of  all  rainstorms  will  create  runoff.  Consequent  problems 
arising  in  the  definition  of  CN  from  field  rainfall-runoff  data  are  discussed. 

Background 

The  Curve  Number  Method  is  a  widely  used  technique  for  estimating  storm 
runoff  depth  from  rainfall  depth.  It  was  pioneered  and  developed  by  the 
USDA-Soil  Conservation  Service  (SCS).  The  primary  source  document  for  the 
method  is  their  National  Engineering  Handbook,  Section  4  "Hydrology",  or  for 
short  "NEH-4"  (14) 

The  basic  equation  is  simply 

Q  =  (P  -  .2S)2/(P  +  .8S)  (P>.2S)  (1) 

Q  =  0  ( P< . 2S ) 

where  Q  and  P  are  the  runoff  and. rainfall  depths  and  S  is  a  linear  index  of 
watershed  storage  transformed  to  the  index  "Curve  Number"  by 

CN  =  1000/(10  +  S)  (S  in  inches)  (2) 

or  CN  =  25400/(254  +  S)  (S  in  mm) 

Curve  numbers  reflect  the  land  condition,  and  tables  and  charts  of  CNs 
defined  according  to  soil  type,  land  use,  cover,  and  watershed  moisture  status 
are  given  in  numerous  SCS  and  other  agency  documents. 

The  watershed  antecedent  moisture  condition  (AMC)  is  one  of  the  most 
influential  factors  in  determining  CN.  NEH-4  gives  in  its  Table  10.1 
conversions  through  three  moisture  classes,  AMC  I  and  AMC  II,  and  AMC  III, 
where  condition  II  is  the  benchmark  situation  upon  which  a  watershed  is 
described,  and  condition  I  and  III  are  the  "dry"  and  "wet"  situations.  An 
abbreviated  form  of  this  relationship  is  given  here  at  Table  1.  NEH-4  also 


^Professor,  Department  of  Forest  Resources,  and  Professor,  Department  of 
Civil  and  Environmental  Engineering,  Utah  State  University,  Logan,  Utah 
84322. 

2Hydraulic  Engineer,  U.S.  Department  of  Agriculture,  Agricultural  Research 
Service,  Watershed  Research  Unit,  Columbia,  Missouri  65203. 

Graduate  Student,  Agricultural  University,  Wageningen,  Netherlands. 
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gives  suggested  approximate  five-day  antecedent  rainfall  depths  necessary  to 
achieve  these  three  categories  (which  are  not  given  here). 

Table  1.  Relationship  Between  CN  and  AMC  Antecedent  Moisture  Class 


I 

II 

III 

100 

100 

100 

87 

95 

98 

78 

90 

96 

70 

85 

94 

63 

80 

91 

57 

75 

88 

51 

70 

85 

45 

65 

82 

40 

60 

78 

35 

55 

74 

31 

50 

70 

22 

40 

60 

15 

30 

50 

9 

20 

37 

4 

10 

22 

-0- 

-0- 

-0- 

Source:  NEH-4,  Table  10.1  (condensed) 

Generalized  Relationships 

Equation  1  describes  a  family  of  storm  runoff  curves  in  the  parameter  S 
(or  its  tranform,  CN).  It  is  easily  put  into  dimensionless  form  by 
standardizing  on  S,  yielding 

q/s  =  (P/S-.2)2/ (P/S+.8)  (P/S>.2)  (3) 

q/s  =  q  (p/s<.2) 

Also,  the  information  on  the  AMC-CN  relationship  Table  1,  may  be 
abbreviated.  The  S  values  corresponding  to  the  CNs  (as  given  in  Table  1 
herein  or  NEH-4  Table  10.1)  are  found  to  present  a  close  correlation: 


Si  =  2.281  Sn 

r2  =  0.999 

Se  =  0.206  in 

(4) 

Siii  =  0*427  Sjj 

r2  =  0.994 

Sp  =  0.088  in 

(5) 

The  above  equations  pertain  to  the  range  50<CNjj<95,  which  encompasses 
most  estimated  and  experienced  CNs.  Because  Condition  II  is  the  reference 
status  in  Eqs.  1  and  3,  substitution  of  2 . 281Sj j  and  0 . 427S j j  (Eqs.  4  and 
5)  into  3  leads  to 


Q/S  =  (P/S-0.085)2/(P/S+0.342) 

P / S> . 085 

(6) 

Q/S  =  (P/S-0.456)2/ (P/S+l .824) 

P / S> . 456 

(7) 

3 


AP-3 


where  (6)  and  (7)  yield  AMC  III  and  AMC  I  runoffs,  and  S  is  defined  on  AMC  II. 
These  three  functions  are  shown  in  Figure  1. 


For  more  clarity,  the  rainfall  runoff  equations  may  be  grouped  and  written  as 


Qin  =  (P-.085SII)2/(P+.342SII) 

(8) 

Qn  =  (P-.2Sn)2/(P+.8Sn) 

(1) 

Qj  =  (P-.456Sn)2/(P+1.824$n) 

(9) 

In  the  above,  the  effective  term  rainfall  (inside  the  parentheses  in  the 
numerator)  must  be  positive.  Q  =  0  otherwise.  Parenthetically,  it  should  be 
noted  here  from  equations  4  and  5  that 


(SI/SII>  s  (SII/SIII}  a  2'28 


(10) 


Also,  direct  algebraic  conversion  of  AMC  is  possible  via  Eqs.  1,  4,  and  5. 
Substitution  and  simplification  leads  to 


CN 


I  = 


CN 


II 


osnrtMTTfT 


II 


r2  =  0.996 

=  1.0  CN 
e 


(ID 


8^II 

CNlll  =  TJT427+n0^57TTN^ 


r2  =  0.994 

=  0.7  CN 
e 


(12) 


Equations  11  and  12  are  similar  to  equations  suggested  previously  by  Sobhani 
(12)  and  are  similarly  derived. 

Probability  Descriptions 

A  broader  interpretation  of  the  AMC  categories  is  as  "error  bands"  or 
envelopes,  indicating  the  experienced  variability  in  rainfall-runoff  relation¬ 
ships.  In  this  role,  site  moisture  per  se  acts  as  a  surrogate  for  all  other 
sources  of  variability,  such  as  rainf all  intensity,  spatial  variability  of 
watershed  or  storm  properties,  or  basic  shortcoming  in  the  structure  and 
coefficients  in  Eq.  1.  This  role  for  AMC  is  found  originally  in  NEH-4,  and 
has  been  more  recently  elaborated  by  Rail i son  (7),  Rail i son  and  Cronshey  (8), 
Rallison  and  Miller  (5),  and  Springer  et  al .  (13). 

Hjelmfelt  (5)  found  that  the  established  AMC  relationships  fairly 
accurately  described  the  90  percent,  50  percent,  and  10  percent  cumulative 
probability  of  the  runoff  depth  for  a  given  rainfall  via  the  CN  method, 
corresponding  to  AMC  I,  II,  and  III,  respectively.  An  illustration  is  given 
in  Figure  2.  Thus,  equations  1  and  6  to  9  have  runoff  probabilities  associated 
with  them.  For  example,  if  at  AMC  I,  there  is  then  a  10  percent  probability  of 
runoff  depth  Q  exceeding  that  indicated  by  equation  9.  If  say,  CNjj  =  80  and 
P  =  2  inches,  runoff  expectations  would  then  be: 

90%:  Q>0 . 11  in 

50%:  Q>0.56  in 

10%:  Q>1 . 12  in 

where  the  above  calculations  were  made  from  Eq.  9,  1,  and  8,  respectively. 
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This  may  be  seen  generally  as  a  conditional  probability  distribution  at  a 
rainfall  P  with  a  mass  on  the  P  axis,  for  which  Q  =  0.  This  probability  mass 
is  the  fraction  of  occurrences  when  the  effective  rainfall  is  zero,  or  when  P 
is  less  than  or  equal  to  the  initial  abstraction. 

This  occurs  when  P<0.085S  for  condition  III,  when  P<0.2S  for  condition  II, 
and  P<0.465S  for  condition  I.  From  this  reasoning,  and  obvious  behavior  at 
the  extremes  of  P  =  0  and  P  =  =>  ,  Table  2  may  be  established. 

A  log  normal  distribution  closely  approximates  Table  2,  with  a  mean  (yn) 
of  -1.609  and  standard  deviation  (a)  of  0.67.  The  transformation  between  AMC 
values  can  then  be  expressed  in  continuous  and  graphical  form  as  shown  in  Fig. 
3.  Given  the  nature  of  the  "data"  from  which  this  was  derived,  it  should  be 
considered  approximate. 


Table  2.  Runoff  Probability  for  P/S* 


p73 

pfrarsan - 

Comments 

0 

0 

0.085 

0.10 

AMC 

III 

0.20 

0.50 

AMC 

II 

0.456 

0.90 

AMC 

I 

00 

■1 

*S  defined  on  AMC  II 


Storm  Size  Perspectives 

There  is  largely  undocumented  undercurrent  of  concern  about  the  storm  size 
requirements  for  application  of  the  CN  method.  A  main  thread  of  the  original 
intent  was  flood  event  design,  and  early  development  was  through  annual  maximum 
daily  rainfall  and  runoff  data.  This  suggests  some  selection  of  larger  events, 
though  the  precise  scale  of  "largeness"  eludes  definition.  Smith  and 
Montgomery  (11)  claim  that  "P/S<0. 4. . .poorly  (if  at  all)  define  a  CN  in  any 
case...",  though  they  gave  no  reasons.  Hjelmfelt  (4)  suggested  that  a  "small" 
storm  is  one  which  "...P<0.2S  where  S  is  determined  from  the  AMC  I  curve 
number,"  and  cautioned  users  "...  to  eliminate  small  events  from  the  data  set." 
As  developed  in  preceding  sections  this  smallness  threshold  suggested  by 
Hjelmfelt  is  P<0.456S,  or  Pr(Q/S>0)  =  0.90,  where  S  is  defined  on  AMC  II.  In 
both  of  the  above  cases,  the  concern  was  determination  of  CN  from  field  data, 
and  in  the  latter  case,  possible  biasing  of  the  smaller  scale  data  due  to 
including  only  the  events  in  which  runoff  occurred. 

From  the  preceding  development,  it  might  be  seen  that  Hjelmfelt's  "small" 
storm  is  actually  "large"  in  terms  of  runoff  probability.  Thus,  the  following 
classification  is  proposed,  defining  relative  storm  size  on  CN-based  runoff 
probabilities. 
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Table  3.  Relative  Storm  Sizes 


P  /  S 1 1 


Prob. (Q/Sji>0) 


Category 


From  To 


From  To 


0  0.085 
0.085  -  0.20 
0.20  -  0.456 
0.456  -  0.60 


0  0.10 
0.10  -  .50 
0.50  -  .90 
0.90  -  .95 


Very  Small 


Smal  1 
Medi urn 
Large 


>0.60 


>0.95 


Very  Large 


Using  this  as  a  guide,  the  user  caution  for  CN  determination  from  data 
could  be  recast  as  "use  only  large  or  very  large  events",  or  more  precisely, 
"use  only  events  for  which  P/Sji>0.46  or  Pr  (Q/Sji>0)>0.90.  This  is  in 
harmony  with  Smith  and  Montgomery's  advice  of  avoiding  P / S j i<0 . 4 . 

This  pursuit  says  nothing  about  the  probability  of  P/Sji,  which  are 
geographically  dependent.  For  example,  for  a  near-impervious  parking  lot  of 
CN  =  99,  S  =  0.10  in,  and  almost  all  storms  would  be  "large"  events.  On  the 
other  hand,  on  a  very  porous  forested  watershed  of  CN  =  40,  S  =  15  inches,  and 
a  storm  of  3  inches  would  be  "small".  The  occurrence  of  P / S  j  j  is  largely  a 
matter  of  geologic  and  meteorologic  happenstance. 

Determining  CN  from  Data 

Although  the  CN  method  is  usually  used  in  a  prediction  mode  with  CN  values 
estimated  by  handbook  methods  from  soils  and  vegetation  information,  it  should 
also  possible  to  determine  CNs  from  observed  event  rainfall  and  runoff  data 
for  small  watersheds.  Equation  1  can  be  solved  for  S  to  yield 


(13) 


Thus  any  P  and  Q  pair  with  Q>0  leads  to  an  S  and  to  a  CN  via  Eq.  2.  From  the 
rhetoric  in  the  preceding  section,  however,  low  P/S  events  should  be  excluded 
when  determing  CN  from  field  data.  The  ideal  sample  should  contain  a  balanced 
random  scatter  of  events  on  both  sides  of  AMC  II,  which  would  happen  only  in 
high  P/Sji  situations.  Low  P/Sn  events  would  have  runoff  only  above  AMC 
II.  Therefore,  analysis  to  determine  a  Sj i  as  a  mean  (or  median)  of  a  series 
of  events,  using  Eq.  13,  should  exclude  all  points  for  which  P/Sji<0.46. 

This  is  awkward,  because  Sn  is  defined  or  estimated  from  the  same  data  body 
used  to  determine  sample  point  acceptability.  The  task  is  to  censor  the  data 
set  so  that  all  remaining  points  used  to  calculate  Sn  (as  a  mean  or  median) 
have  P/Sjj  >0 . 46 . 
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This  may  be  accomplished  by  tri al-and-error  on  a  computer.  The  procedure 
is  briefly  described  as  follows: 

1)  Use  the  biggest  rainstorm  and  calculate  S  from  Eq.  13,  and  CN  from 
Eq.  2. 

2)  Check  for  P/S>0.46. 

3)  If  P / S>0 . 46  then  add  next  biggest  storm  to  the  calculation,  and  use 
mean  S  values.  Go  back  to  step  2. 

4)  Include  all  events  down  to  point  where  the  last  P  divided  by  the  mean 
(or  median)  S  is  greater  than  0.46. 

To  illustrate  this,  a  number  of  small  watershed  data  sets  are  described  in 
Table  4.  While  Table  5  shows  the  outcome  of  the  CN  determinations.  Figure  4 
shows  a  plot  and  fit  by  this  "P/S"  method  for  the  Chickasha,  Oklahoma  data  set. 

The  exercise  leads  to  some  interesting  points.  First,  the  "P/S>0.46"  data 
censoring  produces  lower  CNs,  as  shown  in  Table  5.  This  occurs  because  the 
high  AMC  events  at  low  P  infer  higher  CNs,  and  these  are  excluded.  In  a  larger 
data  treatment  of  90  watersheds  (not  included  here)  this  was  the  finding  in 
every  case. 


Table  4.  Watershed  and  Data  Description 
Name  Location  AreafAC)  Events  PerToJ  Source  oi-§at a 


3  Bar  D 

Ariz. 

82 

171 

1956-79 

J.D.  Hewlett, 

ref.  3 

(USFS  data 

Badger  Wash  2A 

Colo. 

107 

66 

1953-66 

USGS,  via  ref 

.  2 

Dugout  Cr. 

Wyom. 

518 

37 

1965-72 

USGS,  via  ref 

.  10,  and 

J.  Rankl 

1982-83 

Ephiriam  A. 

Utah 

11.2 

67 

1916-29 

USFS,  via  ref 

1 

Chickasha  5 

Okla. 

23.7 

73 

1967-77 

USDA,  ARS 

Zululand  16 

S.  Africa 

i  796 

43 

1976-79 

ref;  6 

Table  5 

Data  Analysis 

and  Curve 

Number  Determination 

All  Events 

P/S  >  .461 

Watershed 

tN 

. Rl nl  . 

N 

•£n  - 

FHn) 

N 

3  Bar  D 

58.5 

2.04 

171 

26.3 

14.14 

1 

Badger  2A 

93.8 

0.44 

66 

90.7 

0.76 

22 

Dugout 

95.0 

0.49 

37 

93.6 

0.65 

18 

Ephriam  A 

89.0 

0.44 

67 

undef^ 

0 

Chickasha  5 

77.6 

1.37 

73 

68.6 

2.87 

9 

Zululand  16 

76.1 

2.07 

43 

65:8 

4.95 

10 

Notes:  1.  The  parameter  S  here  is  the  mean  S  in  all  the  events  included.  It 
is  taken  to  represent  AMC  II. 


2.  No  events  whith  P/S  >  0.46.  Maxium  storm  depth  in  data  set  is  1.50 
inches. 
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Second,  not  all  events  are  useful  in  detemining  CN.  As  an  extreme,  the 
Ephiram  A  data  set  contained  no  events  for  P/S  >  0.46,  and  thus  CN  could  not 
be  defined.  This  problem  is  minimized  on  high  runoff  sites  with  heavy 
rainstorms.  At  the  alternate  extreme,  low  response  sites  such  as  forests  on 
deep  soils  in  areas  of  tranquil  rainstorms  might  require  several  decades  of 
data  to  include  a  high  P/S  event.  As  an  example,  a  site  with  CN j j  =  50  would 
require  a  storm  of  4.6  inches  to  achieve  P/S  =  .46.  In  widespread  areas  of  the 
western  U.S.,  daily  rainfalls  of  4.6  inches  are  at  a  return  period  of  at  least 
several  hundred  years.  Thus,  lower  CNs  may  defy  practical  identification 
within  our  lifetime. 


Third,  while  this  "P/S>.46"  procedure  censors  data  to  minimize  calculation 
bias,  it  also  reduces  effective  sample  size  and  thus  adds  uncertainty,  which 
is  a  function  of  N~0*5.  Some  practical  accomodation  of  these  two  problems 
would  be  useful,  though  none  is  offered  here.  The  use  of  0.46  as  a  rejection 
level  for  P/S  is  based  on  judgment  only  and  draws  from  the  90%  probablility 
associated  with  it. 


Last,  the  criteria  of  P/S  >  .46  leads  to  an  alternate  expression.  By 
manipulation  of  Equation  3,  it  can  be  easily  shown  that 

Q/P  =  (P/S  -  .2)^/ ( (P/S) (P/S  +  .8))  (14) 

Here,  when  P/S  =  .46,  Q/P  =  .12.  This  suggest  that  Q/P  >  .12  might  be  used  a 
convenient  rule-of-thumb  for  quick  checks  of  data  sets.  At  least  ene  point 
should  have  Q  >  .12  P  for  acceptable  (P/S  >  .46)  CN  determination.  This  is  a 
necessary  but  not  sufficient  condition. 
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Figure  Titles 


Figure  1. 

The  CN  method  rainfall -runoff  relationships,  standardized  on  the 
storage  parameters  at  AMC  II.  The  upper  and  lower  curves  are  for 
AMC  III  and  I,  respectively.  The  figure  represents  equations  3,  6, 
and  7. 

Figure  2. 

Experienced  variability  in  CN  and  handbook  AMC  I  and  III  criteria 
from  Hjelmfelt,  Kramer,  and  Burwell  C). 

Figure  3. 

Probability  of  Q/ S jj >0  as  a  function  of  P / S j j .  This  fits  a 
lognormal  distribution  closely. 

Figure  4. 

Storm  rainfall  and  runoff  for  Chickasha  Watershed  5,  Oklahoma. 

Only  the  points  as  "+"  are  used  in  the  calculation  of  CN  =  68.6. 

In  this  case,  the  inclusion  of  the  next  highest  point  (P  =  2.27  in, 
Q  =  0.06  in)  would  lead  to  CN  =  67.4,  S  =  4.84,  and  P/S  =  .44, 
approximately. 
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Appendix  2.  Examples  of  Job  Control  Language  for  the  ORSER  programs 
on  the  IBM  370  computer,  and  their  outputs. 


TPINFO  program: 

****  tso  foreground  hardcopy  **** 

DSNAME=SEANH15. AD.OR.TPINFO.CNTL 


//SEANHAZ 1  JOR  (4071  090299* RJ029) * • ZFVENBEBfiEN «  * MSGLEVEL = (1*1)* 
//  CLASS=B*TIHE= ( *20) *PRTY=3 
/♦ROUTE  PRINT  RMT29 

//JOBLIB  DO  DSN=SFANH02.ATJJ,ORSERLIR*DISP=SHP 
//STEP  1  EXEC  PGM  =  TPINFO*REGION  =  20OK 
//FT06F001  DD  SYSOUT  =  A 

//FT09F001  DD  SYSOUT=A*DCB= (RECFM=FA*RLKSIZE=133) 

//FT08F001  DD  UN  I  T  =  T  APE9  *  VOl.=SER=WOB5P5  *  L  ABFL=  (1 1  ♦  SL  *  *  IN )  * 

//  DCB= (RECFM=VBS*LRECL=3696*BLKSIZF=3700) * 

//  DSN=SEANH02. AT JJ.0R.0K6-74.S4.DAT A *DISP= (OLD* KEEP) 
//FT05F001  DD  * 

ID 

table 

response 

Off  ir  f  fOO  of  mo  T£  SfNflNf.  OF  f*BTM  orsOuwcES 

OutDUt  *  7mF  pCnnSvlvan!  a  S  T  4  T  F  u*|v£Osi  TY  - 


TAOf  InFOBMATJON 


FILE  lOCNTlf ICATION  BfCOPO 


Tape  SffJl>e*CI*G  nu"HC«S. 
«SS  ADJUSTED  L I WE  Lr*GTH 
0AT£  of  E«POSUPE . 

iNAGf  FORMAT  CEmTFM 

LaT I TUOC . 

LONG! TUOf . 

NA0I®  POINT 

L  AT  I TUOC . 

LONGITUDE . 


TAPE  IDENTIFICATION . . .  BS0016  •  FILE  1 

fltgmtlinc  oo  usfb  identification . 

DL  •  T  F  own  Of  SC®  T  PT  I  ON . . . 

SCannCR  DC  SC® I PT I  On . 

OATf  This  OATA  SET  -AS  . . 

USE  p  MHO  POfPAOen  THIS  OATA  sct . 


SUBSET  SOUBCF  TAPC .  .  FILE  1 

SCFnC/FOahc  . .  1706-1633? 

DA»  *iUH»FO  BfLATlVC  TO  LAUNCH 

A  T  MHf  OF  OMSfHV  A  T  I  ON . .  706 

HOOP  at  T  I  Hf  OF  OHSEQVATION . . .  16 

-INUTC  AT  T|Nf  OF  opSCBvaTIOn .  33 

TfMS  of  «CCON0S  AT  T|«C  of  ORSfOVATlON .  ? 


A  OF  A 
3240 

?N  .HIM*  74 

N  -  34  0C0  40  "IN 

W  -  OQP  OF  G  31  H I N 

N  -  3*  flFO  *ft  -I" 

W  -  OOH  OCG  1®  "IN 


OATA  BECOPO  length.. . 

SUN  FLEVATION  angle . 

Sun  azimuth . . 

SATELLITE  HFA01MG . 

OfvOLUTTON  NUmREB.... . 

TtP£  OF  OHHIT . 

COTS  MISSION  nU"BCO . 

niorcT  no  occooneo  mss  data 

HSS  OATA  ACQUISITION  SITF.. 


channel  NUNMfOS  AMO  SPFCTOal  Bands 


Channel  Nii-Pfo  SPECTRAL  aanO  I  m  I  C°OnS  ) 

1  0.S0  -  0.60 

2  0.60  -  0.70 

.1  0.70-0.60 

4  0. AO  -  I  .  10 


32®6 
60  nfG 
I  0 a  OCG 
I  69  DCG 

4644 

n 

l 

D 

G 


TamlE  OF  CONTENTS  BfCOOD 


mfg.  LInC 
I 


Cno  l Inf 

l  »  oo 


Hf  0  •  EL- 

?*?S 


End  el- 

?HOO 


t  INC  INC. 


SCAN  LINES  CLEMENTS 
1100  376 


FI»#ST  MfCOUfi  IN  This  fjlE: 


AP-11 


SUBSET  program: 

****  tso  foreground  hardcopy 
DSNAMF=SEANH15.AD,0R.SUBSET.CNTL 

//SEANHAZ1  JOB  (4071090299tRJ029) t • ZE VENRERGEN » ♦ MSGLE VEL= ( 1  *  1 ) ♦ 
//  CLASS=B»TIME= ( *20) *PRTY=3 
/•ROUTE  PRINT  RMT29 

//J0BLI8  DO  DSN=SEANH02.ATJJ.0RSERLIRtDISP=SHR 
//STEP  1  EXEC  PGM=SURSET,PEGION  =  300K 
//FT06F001  DD  SYS0UT=A 

//FTOflFO 0 1  DD  UNIT  =  TAPE9» V0L  =  SFR  =  W085fl5» DISP= (OLD*KEEP)  * 

//  LABEL=  ( 1 1  *SL  *  *  IN)  * DSN  =  SE ANHO 2.  ATvJJ. OR. 0K6-74.S4.DAT A 
//FT09F001  DD  UN  I T  =  T APE9 * D I SP= ( NEW , KFFP , DELETE )  » 

//  DSNsSE ANH 1 5. SUBSET . OKLAHOMA. DAT A* 

//  DCB= (RECFM=VBS*LRECL=3696,0LKSIZE=37OO) .LABFL= ( 1 tSL) 
//FT05F001  DD  * 

BLOCK  400  500  2501  2600  1  1 

CHANNELS  1234 

END 

n  OFFICE  FOP  wEmOTE  SENSING  OE  EARTH  RESOURCES 

Output:  THE  PENNSYLVANIA  STATE  UNIVERSITY 


SUBSET  PHOGHam 

CSIZ«C>r3SSSll3CSXC3S3K?3S* 

subset  erom  TAPE:  .  file  1 

SUBSET  ONTO  TAPE:  USOOie.  ERE  1 

CHANNELS  ON  TAPE:  1  2  3  A 

CHANNELS  SUBSET!  1234 


INPUT  Tape  TABLE  OF  CONTENTS 


BEG.  LINE 

END  LINE 

BEG. EL. 

ENO  EL. 

LINE 

INC. 

EL.  INC. 

SCAN  LINES 

ELEMENTS 

1 

lion 

2425 

2600 

1 

1 

1 

1100 

376 

reouesteo  table 

OE  CONTENTS 

BEG.  LINE 

ENO  LINE 

BEG. EL. 

ENO  EL. 

LINE 

INC. 

EL.  INC. 

Scan  LINES 

ELEMENTS 

400 

son 

2501 

2IS00 

I 

1 

0 

0 

OUTPUT 

tape  table 

OE  CONTENTS 

BEG.  LINE 

ENO  L INE 

0FG.EL. 

EnO  EL. 

LINE 

INC. 

EL.  INC. 

SCAN  LINES 

ELEMENTS 

400 

son 

2501 

2f.no 

1 

I 

1  n  1 

1  00 

NMAP  program: 


****  ISO  FOREGROUND  hardcopy 
DSNAHE=SEANH15.AD.0R.NMAP.CNTL 

//SEANHA21  JOB  (4071090299, PJ029) ♦ • ZFVENBERGEN • *MSGLEVEL= ( 1 , 1 ) 
//  CLASS=B*TIME= ( *20) ,PRTY=3 
/*R0UTE  PRINT  PMT29 

//JOBLIB  DD  DSN=SFANH02. AT J J . ORSERL I R ♦ 0 1 SP  =  SHR 
//STEP  1  EXEC  PGM=NMAP» REGION  =  500K 
//FT06F001  DD  SYSOUT=A 

//FT08F001  DD  UNIT=TAPE9» VOL=SER=W06B24, 

//  0SN=SEANH15. SURGM. OKLAHOMA. DATA, DISP= (OLDf KEEP) ♦ 

//  DCB= (RECFM=V8S,LRECL=3696,BLKSIZF=3700) ,LABEL= ( 1 tSL, • IN) 
//FT09F001  DD  S YSOUT= A ♦ DCR= < RECFM=F A , RLKS I ZE= 1 33 ) 

//FT05F001  DD  * 

BLOCK  401  450  2S01  2600  1  1 

channels  2  4 

end 


Output  NMAP  program 


* 


OFT  ICE  FOB  REMOTE  SENSING  OF  EARTH  RESOURCES 

The  Pennsylvania  state  university 
AT  01  0 


TAPE  name .  RS001B 

SCENE/FPAME  ID .  1706-16332  FLIGHTLINE  OR  USER  ID . 

DATE  OF  EXPOSURE .  29  JUN,  74  TAPE  SEQUENCING  NUMBERS . 4  OF  4 


BRIGHTNESS 


CLASSIFICA 


T  I  0  N  “  A  P  P 

33X::=c:sri:::::3i 


I  N  G 


CHANNELS  USED:  2  4 

NO.  OF  GRAY  SCALE  LEVELS  FOR  EACH  CHANNEL! 
CHANNELS:  2  4 

120.  120. 

CLASS  LIMITS  DETERMINED  AUTOMATICALLY. 
CLASS  SPECIFICATIONS 


SYMBOL  LIMIT 


X  21.33 

♦  22.50 

23.01 

100.00 


PROGRAM:  nmap 

OATE  AND  TIME:  0:  0 

TAPE  name:  RSOOlfl 

BLOCK  SPECIFICATIONS!  401  450  2501  2600  1 


SYMBOL 

COUNT 

SYMBOL 

LIMIT 

COUNT 

PER  CENT 

X 

21.33 

1205. 

24. 

♦ 

22.5B 

1320. 

26. 

- 

23.01 

1303. 

26. 

100.00 

1172. 

23. 

average  response  per  ELEMENT  FOR  EACH  CHANNEL: 
CHANNELS:  2  4 

33.0  23.7 


MAXIMUM  RESPONSE  IN  THE  BLOCK  FOR  EACH  CHANNEL: 
CHANNELS:  2  * 

50.  32. 


MINIMUM  RESPONSE  in  the  BLOCK  FOR  EACH  CHANNEL: 
Channels:  2  4 

18.  12. 

AVERAGE  NORM  per  Element:  40.7 

MAXIMUM  NORM  IN  THE  BLOCK!  50. 

MINIMUM  NORM  IN  the  BLOCK:  30. 
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PROGR 

am  : 

NMAP 

DATE 

AN  0 

time: 

T  aPE 

name 

:  RS001R 

BLOCK 

SPE 

C  I  E  I C  a  T I onS 

8E 

G  INM 

ING  LINE  * 

ENOING  LINF  *50 


BEGINNING  ELEMENT  ?5ni 

ending  element  2600 

LINE  INCREMENT  1 
ELEMENT  INCREMENT  1 


I 
T 
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?S01  1 25 06  I  251  1  1  25 16  12521  12626  125  31 
I  I  I  I  I  I 

*01  I--**-  ~  *XX  —  .XX - — 
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*0*  I**- —  -  -.XX.XXX.-  -» ♦  — X 

*05  I  —  —  -.xxxxx. - *  —  XX- 

*06  I  -  -xxxx.xx .xxxxxx 

*0?  I  -  -...X*.-  -.x.xxx 

*08  I  -  — ...X»...~. .XXXXX 

*09  I  -*-  —  -.xxxxx. — x.xxxx 

*10  I  -  — ..XXX— ..XXXXX— X 

*11  I  *♦.-  — .-.xxxx.xxxxxxx- 

*12  I  — ♦  ♦  ♦♦♦xxxxxxxxxx - 

*13  I  ♦  *-  -  -  - .xxxxxxx.-  - 

414  I  —  - -  -  -. X.XXX XXXXX..  - 

*15  I  -  - -  — ♦-♦-xxxxxxxxxx-* 

*16  I  ♦ — .xx*. —  -».*.xxxxxxxx*-x 

*17  I  .-..X.XX - - ..XXXXXX - XX 

*18  I  -  .♦  XXX**  —  —  —  XXXXXXX  -XXX 

*19  I  -  — - — *.*..XXX*-*XXXX 

*20  I  -  -XX* - x*xxxxxxxxxx* 

*21  I  -  -X*XX****— -XX ♦  XXX XX *-*x * 

*22  I  xxxxxx*.-  —  «xxxxx-*— ♦** 

*23  I-  —  - ♦ X  X  X  *  * - **xxxxxxx-  -*♦ 
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*29  I— **\X-  ♦  ♦  — xxxxxxxxxxx— 

*30  I- — x*xx*— xx***xxxxxx*xxxxxx*- 
*31  I - X*— X*  —  *XXX  — **X*X*  — 

*32  I** - *x*-xx— *xxxx*  -x*— ♦— 
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Appendix  3.  Examples  of  SAS  programs  used  in  the  Landsat  analysis. 


****  TSO  FOREGROUND  HARDCOPY  ••** 

DSNAME=SEANHl5.AD.SAS.AVREFL.CNTL 

//SEANHAZ1  JOB  (A071090299.PJ029)  .  *ZEVENRERGEN» ,CLASS=C» 

//  TIME=<0.20) »PRTY=13 

/•ROUTE  PRINT  RMT29 
//  EXEC  SAS 

//OKLAHOMA  DO  UNIT=SYSDA.DSN=SEANH03. AD . D I G I T 3 . 0 AT A . D I SP=SHR 
//SYSIN  DO  * 

Data  oklaj 

Infile  Oklahoma; 

INPUT  DATE  1-3  WTRSHED  $  5-6  PIXNUM  H-9  DIG*  11-12  DIG5  1*-15 
DIG6  17-18  DIG7  20-2 1 l 
IF  DATE=  96  THEN  SIN=0.78t 
IF  DATE  =  2*0  THEN  SIN=0.81I 
IF  DATE=3*8  THEN  SIN=0.**l 
IF  DAT E  =  1 8 1  THEN  SIN  =  0.R71 
IF  0ATE=2B9  THEN  SIN=0.63I 
REFLAs (DIG4/289) /SINI 
REFL5=  (DIG5/306) /SIN! 

REFL6= (0IG6/28*) /SIN; 

REFL7= (0IG7/125) /SIN; 

OUTPUT  ; 

PROC  SORT  oata=okla; 

by  wtrshed  oate; 
proc  means  noprint; 

by  WTRSHED  DATE; 

VaR  REFL*  REFL5  REFL6  REFL7 l 

OUTPUT  OUT=REFMFAN  MEANaRMEAN*  RMFAN5 

RMEAN6  RMEAN7  STD=SO*  SD5  S06  S07J 
PROC  PRINT; 

title  REFLECTIONS  by  WATERSHED  AND  date; 

PROC  PLOT  uniform; 

PLOT  RMEAN4*DATE=*#*  RMFAN5*OATE=»« »  RME AN6*0 A TE = • * • 

RmEan7«OATE=,*,/hZFRO  VAXIS=0.10  TO  0.2B  RY  0.01  HPOS=80; 
Calculation  of  the  average  reflection  in  four  bands  on  five  different 
dates  for  Landsat-1  data. 

•*«•  TSO  FOREGROUND  HARDCOPY  *«•• 
DSNAME=SEANHl5.AD.SAS.REFCHAR.LNDST2.CNTL 

//SEANHAZ1  JOB  (4071090299»RJ029) » •ZFVENBERGEN' «CLASS=C. 

//  TlME=  (0.20) .PPTY  =  3 
/•ROUTE  PRINT  RMT29 
//  EXEC  SAS 

//DIG  DO  UNIT=SYSDA.DSN=SEANH15.DIG.WFDC1 .DATAtDISPsSHR 
//SYSIN  DD  * 

data  dig; 
infile  dig; 

INPUT  LINE  EL  DIG*  DIG5  DIG6  DIG7I 
SIN.0.8481 

REFL4=0.1775» ( (DIG4/63.5) *0.10) /SIN; 

REFL5=0.2074*{ (DIGS/85.81) ♦0.07)/STn; 

REFL6=0.25A0* ( (DIG6/95.49) ♦0.07)/SIN; 

REFL7»0. 1261* ( (D IG7/ 15.71 ) *0. 14) /SIN ; 

DROP  LINE  EL  SIN; 

proc  means; 

VAR  REFL*  REFL5  REFL6  REFL7 l 

title  means  and  variances  of  reflection; 
proc  corr  oatasdig  out=correl; 

VAR  REFLA  REFL 5  RFFL6  REFL 7 ! 

title  correlations  and  means  of  the  four  bands; 

PROC  CHART  DATAsDIG; 

HRAR  REFLA  REFL5  REFL6  REFL7  * 

TITLE  REFLECTION  distribution  in  FOUR  bands; 

Program  prints  horizontal  histograms  for  the  distribution  of  reflection 
values  in  four  bands,  for  Landsat-2  data. 


Calculation  of  reflectance  index  models  for  five  different  dates, 
and  correlation  coefficients  of  these  values  and  ground  truth. 

•***  TSO  FOREGROUND  HARDCOPY  **** 

0SNAHE=SEANH15. ad.sas. VI reel .CORR.CNTL 

//SEANHAZ1  JOB  (4071090299.RJ029) ♦ 'ZFVENRERGEN* .CIASS=C» 

//  TIME=<0.20) ,PRTY=13 

/•ROUTE  PRINT  RMT29 
//  EXEC  SAS 

//OKLAHOMA  DO  UNIT*SY5DA.OSN=SEANH03.AD.OKDIGIT?.DATA.DISP=SHR 
//GROUND  DD  UNlT«SYSDA*nSNrSEANH03. AO. OK, GROUND. DATA»DISP=SHR 
//SYSIN  DD  * 

OATA  okla; 

InfilE  Oklahoma; 

INPUT  date  1-2  WTRSHED  S  4-5  PIXNUM  7-8  OlGA  10-11  DIG5  13-1* 

DIG6  16-17  DIG7  19-201 
IF  DATE=04  THFN  SIN=0.78J 
IF  DATE=0fl  THEN  SIN=0.81! 

IF  DATE  =  12  THEN  5  I NaO ,44  1 
IF  DATE=06  THEN  SIN=0.87! 

IF  DATE=10  THEN  DELETE; 

Al  =  (0IG4/289) /SINI 
A2= (0IG5/306) /SIN; 

A33 (0IG6/284) /STNl 
A4= (0IG7/125) /SIN; 

N06s ( A3-A2) / ( A3*A2) ; 

ND7= ( A4-A2) / ( A  4 ♦ A  2 )  1 
TV  1 6  =  5QRT (ND6* .5)  ; 

T  V 1 7  =  S0RT (N07* .5) ; 

SB  I  =  0.332*A 1  +  0. 6 03*42*0.675* A3* 0.262* A4! 

GV I =-0. 283* Al-0. 66 0*A2*0. 577*43*0.388*44! 

YV I =-0.899* A  1*0. 482*42*0.076* 43-0.04  1*44! 

NS  I =-0.0 16* A  1*0. 131*A2-0.452*A3-0.8R2*A4! 

PV 1 6=SQRT ( ( -2. 50  7-0. 457* A  1 ♦0.498* A3) **2* ( 2. 734* 0 . 498* A2-0 . 543* A 3 > **2> 
PV 1 7  =  SORT ( ( .355*44-. 1 49*42) **2* ( . 355* A2- . 852* A4 ) **2 )  ; 

0VI=2.4*A4-A2; 

AVI=2*A4-a2; 

R5*=A2/A1 ; 

R65=A3/A2; 

R75=A4/A2; 

DROP  0 1 G4  DIG5  DIG6  DIG7  PIXNIJM  SIN; 

PROC  SORT! 

0Y  WTRSHED  date; 

PROC  MEANS  NOPRINT! 

8Y  WTRSHED  DATE! 

VAR  ND6  ND7  T  V 1 6  T V  T  7  SBI  GV I  YVI  NS  I  PVI6  PVI7  DVI  AVI  95* 

R65  R75  A 1  A2  A3  A4! 

OUTPUT  OUT  =  V I ME  AN  MEAN=AVND6  A VN07  AVTVIA  AVTVI7  AVSBI  AVGVI 

AVYVI  AVNSI  A VP V  I  6  AVPV I  7  AVOVI  AVAVI 
AVR54  AVR65  AVR75  AVA1  AVA2  A VA3  AVA4 ! 

DAT*A  ground; 
infile  ground; 

INPUT  rain  1-4  PERCVEG  6-7  FRRIO  9-13  DRBTO  15-19! 

data  correl; 
merge  vimean  ground; 

PROC  PRINT! 

TITLE  NFW  data  SET  CONTAINING  AVERAGE  VEGETATION  indices; 

T I TLE2  and  GROUND  TRUTH  BY  WATERSHED  AND  DATE! 

PROC  CORR  BEST=5I 

VAR  RAIN  PERCVEG  FRRIO  DRRIO! 

WITH  AVND6  AVN07  AVTVT6  A V  T  V 1 7  A VSR  T  AVGVI  AVYVI  AVNSI  AVPVI6 

AVPV I  7  AVOVI  AVAVI  AVR54  AVR65  AVR75  AVA1  AVA2  AVA3  AVA4 ! 


Appendix  4.  Summary  of  Curve  Number  Computations 
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Appendix  5.  Examples  of  some  SAS  and  FORTRAN  programs  used  in  the 
curve  number  analysis. 

Calculation  of  the  mean  curve  number,  and  by  non-linear  regression 
for  natural  data. 

*♦•*  tso  foreground  hardcopy  ♦  •*• 

DSNAMF=SEANh15. AD.SAS.LSQCN.CNTl 

//SEANHAZ1  JOB  (40710902<3R,RJ02R)  *  •ZF'/ENHEBGEN»  ,CLASS  =  Ct 
//  TIMF= (0,20)  , PRT  Y  =  1 3 
/•ROUTE  PRINT  RM T 2B 
//  E*EC  SAS 

//IN  DD  UNIT =SYSDA,DSn=SEANh15,w5. DATA, DISP=ShR 

//STAT  DD  IJNI T=SYSDA,DSN=SEANH1S. STATON. R5. DATA, DISP= (NEW, CATLG) , 

//  SPACE=(TBK.(20»5)»RISE>«OCH=(LRECL=80,RLKSI7E=4S60.RECFm=FR) 
//SySIN  DD  • 

DaTa  in i 
InFTLE  In* 

INPUT  P  Q» 

S  =  5« (P* (2*0) - ( (A* <Qo«2> ) ♦ (5*P*0) ) **0.5) * 

RCN=1  000/ ( 10*S>  t 
OROP  S* 

PROC  PRINT* 

PROC  MEANS  N  MEAN  STD  SUM* 

OUTPUT  OUTaSTATS 

N=NUM  MEANrMEANP  MEANO  MFANCN  STnsSTDP  STOO  STOCN 
SUMsSUMP  SUMO  SUMCN* 

Data  statj 

set  stats; 

FILE  STAT* 

PUT  NUM  MEANP  MEANO  MEANfN  STOP  STDO  STDCN  SUMP  SUMO  SUMON* 
PROC  NLIN  DATA=IN* 

PARMS  CN=50  TO  100  BY  5* 

P 0=200/CN-2t 
IF  P<= PO  then  do; 

MODEL  0=0* 

der.cn=o; 

end; 

ELSE  00* 

MODEL  Q  = (P-200/CN>2) • (P-200/CN*?) / ( P *80 0/CN-R ) ; 
0ER.CN=(400*(P-200/CN*2)*(P*800/CN-B)/(CN*CN)- 

(P-200/CN*?) * (P-?00/CN*2) ♦ (-800/ (CN*CN) ) ) / 

( ( P* BO O/CN-8 ) • (P*800/CN-fl) ) * 

End; 

OUTPUT  OUT  =NEW  PRFDICTEOsQP  RESIDUALsPP* 

PROC  PRINT* 

PROC  'PLOT; 

PLOT  Q*P  QP«P=»*»  /  OVERLAY  HPOS=B0  VPOS=40* 


Calculation 
for  ordered 


of  the  mean  curve  number, 
data. 


and  by  non-linear 


regression 


****  tso  foreground  hardcopy  *•** 

DSNAME=SEANHl5. AD. S AS. LSOORO.CN TL 


//SEANHAZ1  JOB  <4071090299, PJ029) , 'ZEVENRERGEN* ,CLASS=C. 

//  time =(o,2o) , ppt y= i 3 

/•ROUTE  PRINT  RMT29 
//  EXEC  SAS 

//IN  DO  UNI T=SYSDA,DSN=SEANHI5. R5. DATA, DISP=SHR 

//STAT  DD  UNIT  =  SYSDA, DSN  =  SEANH15.MEAN0P0.R5. DATA ,DISP= (NEW, CATLG)  , 
//  SPACE=ITRK,(2O,5).RLSE),DCB=(LRECL=HO,0LKSIZE=456O,RECFm=FB) 
//NEW  DD  UNIT=SYSDA,DSN=SEAMH] 5. LSCN0PD.R5. DATA, DISP= (NEW, CATLG) , 
//  SPACE=(TPK*  (20.5)  .Rt.SF)  ,  0CR  =  ( LRECL  =  8  0  ,  RLKS  I  ZF=4S60  .  RECFm  =  FB  ) 
//SySIN  DO  • 

Data  in) 

Infile  in; 

Input  p  o; 
pboc  sort; 

by  p; 

Data  a ap ; 

SET  IN) 

DROP  p; 

PROC  SORT) 

BY  Qt 

Data  noot; 

SET  ini 
OROP  0) 

Data  MIES) 

MERGE  NOOT  aap; 

S*5* (P* (2*0) -( (4*(Q**2) ) ♦ (5*p*0) )**0.5) » 

RCN=1000/  (  1 0 ♦S )  I 
DROP  SI 
PROC  PR  TNT  I 

PROC  MEANS  N  MEAN  STD  SUM! 

OUTPUT  OUTaSTATS 

NaNUM  ME  ANaME  ANP  MF  AND  MEANCN  STOaSTOP  STDO  STDCM 
SUMaSUMP  SUMO  SUMCNI 

OaTa  stati 

SET  STATS! 

FILE  STAT! 

PUT  NUM  ME  ANP  MEANO  MEANCN  STOP  STDO  STDCN  SUMP  SUMQ  SUMCNI 
PROC  NLIN  DATAsMIES) 

PARMS  CN=20  TO  100  BY  5! 

*  P0  =  200/CN-2I 
IF  P<=P0  THEN  DO! 

MOOEL  0*01 


DER.CNaOl 

END! 

ELSE  DO! 

MODEL  0= (P-200/CN*2) • < P-200 /CN*2) / (P*800/CN-R) l 
DER.CN=(400*(P-200/CN*2)»(P*BO0/CN-fl)/(CN*CN)- 

(P-200/CN+2) • (P-?00/CN*2) • (-HO0/ (CN*CN) ) ) / 

( (p+Bon/CN-tt) * ( p  *  8  o  o /Cn— b ) ) ; 

end; 

OUTPUT  OUTSOLD  PREDICTEO=0P  RESIDHALsHP  PAPMSaLSCN! 
Data  NEW) 

SET  OLD) 
file  new; 
put  p  o  lscn; 

PROC  PRINT) 

PROC  PLOT) 

plot  o«p  qp<.p=»*.  /  overlay  hpos=ro  vpos*4o; 
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Calculation  of  curve  numbers  using  only  the  events  that  meet  the  selection 
criterion:  P/S  >0.46  . 

♦••*  tso  foreground  hardcopy  •*** 

DSNAME*SEANHl5. AD.CN.CNTL 

//SEANHAZ1  JOR  (4071090299. RJ029) ♦ •ZF^ENBERGEN* ,CLASS=C. 

//  TIME* (0*20) ,PRTY=13 
/•ROUTE  PRINT  RmT 29 
//STEP  1  EXEC  FORTGCLG 
//FORT . SYS  I N  DD  * 

DIMENSION  P(37)*Q(37).CN(37)«S(37)*FACT(37) 

integer  N 

HEAD  (8.10)  (P ( I )  ,0  ( I )  *  1*1 *37) 

10  FORMAT  (F4.2.FS.2) 

J  =  37 

WRITE  (6*15) 

15  FORMAT  ( lHl  »///.20X. »P ( I )  *  *  1  OX.  iQ (  I )  •  * 10X. *S ( I )  •  » 1  OX* 

1 »CN ( I )  ♦  ,//) 

DO  100  1=1. J 

S ( I ) =5.* (P ( I ) *2.*0  < I ) - (4.*0  < I ) *0  < I ) *5.*P ( T>  *Q  < I ) >  **0.5> 

CN ( I ) =1000./ (S ( I ) *10. ) 

WRITE  (6.20)  (P(I)  .0(1) *S(I) . CM ( I ) > 

20  FORMAT  (20X.FA.2. 10X.F4.2.09X.F5.2* 10X.F5.2) 

100  CONTINUE 
150  CONTINUE 
ST0T=0. 

DO  200  1=1. J 
STOT  =  ST 0 T ♦ S  (  I  ) 

200  CONTINUE 

SM£ AN=ST0T/J 
DO  400  1=1 .J 
PACT  < I ) =P ( I ) /smean 
IF  (FACT ( I > .LT.0.46)  GO  TO  300 
400  CONTINUE 
GO  TO  500 
300  J=J-1 

GO  TO  150 
500  CONTINUE 

WRITF  (6.15) 

WRITE  ( 6  »  ?0 )  (P ( I ) .0 ( I ) *S ( I ) *CN ( I ) *  1  =  1 . J) 

CNSUM=0. 

DO  600  1=1. J 
CNSUM=CN ( I ) *CNSUM 
600  CONTINUE 

CNMEAN=CNSUM/j 

WRITE  (6.30)  CNMEAN 

30  FORMAT  (//////. 20X. 'THE  AVERAGF  CN  IS  »»F5.2) 

STOP 

END 

//GO .FT08F001  DO  UN T T =S YSO A « DSN=SE ANH 1 5. AD.OUGOUT .UATA«0ISP=SHR 
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Appendix  6.  Location  of  Tombstone  W4  on  line  printer  map, 


program:  nmap 
date  AND  TIME: 
Tape  name:  rsooip 


o :  0 


BLOCK  SPECIFICATIONS 
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ENDING  LINE 
BEGINNING  ELEMENT  2901 
ENOING  ELEMtNT  3000 
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I-4X4XXXXXXX44444 - 44X444XXXX44X 

I  44XX  XX  XX  44  4X4  4  4—  4  4  4  4  4  X  X  X  X  X  X  4  4  X  X 
IXX44XXX444-XXX44XXX-4XXXXX44444X 

I  XXXX  —  44XXXX4XXXXXXXXX4  44  4XXXX4  41 


♦  ♦  X  X  X  4 - — 

44— —44— —4  — 

444444444  —  4  —  —  —  - 

♦♦♦♦♦♦♦XX*  —  « 

--4444-444- 
-4—  X444444 
4  — 4-4XXXX444 
44-4  — 4XXXXXX4 
4  44444  XXX  4-4-  — 

-X*X*XX**X*  — - — 

4  444—  —  4X  XX  XXX  4  4-- - — 

4 XXXX XXXX-* XXX XX 4- — 4-4 — - 


—  44  < 


——44—4— 


IXX4444-4— XXX444X4XXXXXX4*XXXX4-444>—  /  -X4  — XXXXX—  4  X  ♦ 

I  4 4——  —  — 4 XX X X—  X 4 4  4 4 X X  4 4 XX 4XX4 4 4 4 4 4 4 X ♦ 44  4 X X —  444XX44  —  444  —  4  —  —  — 

I - *XXX**X* - XX**XX*4444— 4XXXXX— 4XXX*  44XXXXX44X—  -444-  - 

IXXX44X44XX4XXX444XX4XXX444-44XX - 4X4-XXXX44  -4444X4  - 4  X  X -  - 

!♦♦  —  44XXX*XXX44XX*4XX-44—  4  4  -  4  4  4  X  X  X  X  4  4  X  X  X -4  4  4  4  4  X—  4X4-4-XX4-  -  - 

IXXX-444XXX44XXXX4-444X4-  —  44XX  XXX  XXX  X—  -XX— -XX-  4X44444—44—— 

I4XX44XXXX44XXXXX4XXX  4  — 4  X  — —  ♦ X4  X  X  X  4  4  4 -  -  -  -XXXXXXXXX* - 

I  4  X  4  4XX  4  4  XXXX— —4  X——  4— —  4—  4X4-  -444XX  *-  — ♦  -X  4  X  XX  *4  44  4— 4  4  4- 

IX*XXXXXXX*X,XXX-*XX**XXX*4  —  4X4-XXX4-.4  «  X-  —  44-XXXXXXXXXX—  —  4- 

I4XXXX44444XXX44XXXX  —  4 X X X X X X X X X X X X 4 — —  -4  -  — X  4  X X XX  X 4  4 x X X  X  X X X  4— — * 

IXXX44X-4-XX4-4XXX — —  XXXXXXXXXX1XXX  4—  — — X  444  — 44XXXXXX  — 4XXXXXXXXXX4—  -4 

IXXX44444XX4— XX*X XXXX44X4X444-  -  - X X X - 4 X 4 - X X 4 4 X X 4 4 4 X X X X X X X X X X 4  -4X4 

I-  — — 4  — XXXXXX  — XX  -  — XXXXX44XXX—  -XXXXX44XX4XXX4X4XXX-  -4  X  X  X  *  — 
I  4X44X44  —  XXX4XX44XXX  XXXX  XX4-— XXXX  4  4X4  4XXXXXXX4-  -4  4444XXX4- 

1  4*  —  XXxXXXX**XX*XXXXXXXXXX*-  -XXXX44444XXXXX44  — XX  4  XX  X  X  4  —  4  — 

XX  - —  4  X  X  X  4  4 XX  X  X  4—  — ♦  —  XXXXX44  — 

—  44-4XXXXX—  44XX-  -  ——4  X  X  XX— —  X  4  —  4  4  — 
4  4  X4  4  X  XXXX  XX  4  —  4XX4X4  -  4  X  X  X  ♦  X  4  —  X  X  4  4  4  4- 

XXX  X  4 

xxxxx 

4  4  X  X  4 

I  I  I  I  I  I  I  I  I  I 

29711 


914 

I 

—  ♦ 

-XXXXXX 

915 

I  — — — — ♦ XXX 

XXX  — 

♦XXXX X* 

916 

IX4XX4  — X 

44XXXX4X4XXX 

917 

I  X  X  X  X  4— —4  4 

X  4  4  4  4 

♦  X ♦  — ♦  —  X 

918 

IXXXXXXXXX 

X*- 

♦  X 

919 

IXXXXXXXXX 

XX  XX  * 

-.-♦XX 

920 

I 4  4  xxxxxxx 

X  X  ♦  ♦  4 

-♦XX ♦xx 

I  I 

I 

I  T 

290 1 I 2906 I 29 1 

11291 

6 I 292 1 I 

£  X— — 4  — 4  X  XXXX  — 44  — 4XXX*  - XXX4X4 - 4  — 

(*— — 4XXXXXXXX4XXX4  4-4X44444—  4X4) 

>44X444-UXXXXXXXXX4  -4  XXX  X  X  4—  - — XXX) 
4  X— 4  XXXXX  XXX  4—  — 4XXXX4444  X  X  X  X  4 

I  I  I  I  I  I  I  I 


4—  4X4 
—  4  X  X  4  4  4 
44.4X4  — 
XX—  44- 

I 

29761298 


I  T 
1 29861 2 


AP-24 


Appendix  7.  Parameters  of  linear  and  exponential  relationships  between  curve 
numbers  and  reflectance  index  models  for  9  watersheds. 


RIM 

al 

linear 

bl 

r 

a2 

exponential 

b2 

r 

CH5 

43.82 

246.13 

0.967 

3.894 

3.267 

0.976 

R65 

115.45 

-23.89 

-0.927 

4.849 

-0.320 

-0.944 

ND6 

96.08 

-89.24 

-0.950 

4.588 

-1 .189 

-0.962 

TVI6 

206.54 

-153.84 

-0.955 

6.057 

-2.046 

-0.965 

SB  I 

15.28 

191 .48 

0.938 

3.511 

2.552 

0.950 

i 


